Packages

library(dplyr)
library(ggcorrplot)
library(ggpubr)
library(patchwork)
# library("here")
# library("bookdown")
# library("downloadthis")
library(vegan)
library(plyr)
library(e1071)
library(tidyverse)
library(viridis)
library(GGally)
library(ggrepel)
library(ggordiplots)
library(readr)
library(RColorBrewer)
library(oce) 
library(plotly)
package 㤼㸱plotly㤼㸲 was built under R version 4.0.5Registered S3 method overwritten by 'htmlwidgets':
  method           from         
  print.htmlwidget tools:rstudio

Attaching package: 㤼㸱plotly㤼㸲

The following objects are masked from 㤼㸱package:plyr㤼㸲:

    arrange, mutate, rename, summarise

The following object is masked from 㤼㸱package:ggplot2㤼㸲:

    last_plot

The following object is masked from 㤼㸱package:stats㤼㸲:

    filter

The following object is masked from 㤼㸱package:graphics㤼㸲:

    layout
source("gg_ordisurf_viridis.R")

Data

env_sort <- read.csv(file.path(dataPath,"inputs/LDensShlw_env_sort_2022-09-23.csv")) %>% as.data.frame
otu_sort <- read.csv(file.path(dataPath,"inputs/LDensShlw_otu_sort_2022-09-23.csv")) %>% as.data.frame


table(env_sort$SampID2 == otu_sort$SampID)
< table of extent 0 >
#must be same length and equal to unique sample number length - the ordered lists are assumed to be directly relatable on a row by row basis

dim(env_sort)
[1] 1137  142
dim(otu_sort)
[1] 1137  351

remove any variables with not enough coverage


env_sort <- env_sort %>% select(-c("BO22_lightbotmean_bdmean",
                               "BO22_lightbotltmax_bdmean",
                               "BO22_lightbotltmin_bdmean",
                               "BO22_lightbotrange_bdmean"))

Plot map to check location of samples

env_sort_locations<- ggplot(data = env_sort,
                      aes(x = X.y,
                          y = Y)) +
  theme_classic() +
  geom_point(aes(colour = bathy),
             size = 2) +
  scale_colour_gradient2(low = "red",
                         mid = "yellow",
                         high = "green") +
  ggtitle("Location of samples in sorted env file")

env_sort_locations

Data cleaning

## Removing NAs
otuCompl <- otu_sort[complete.cases(env_sort[, -c(1:23)]), ]
envCompl <- env_sort[complete.cases(env_sort[, -c(1:23)]), ]

## Removing observations with less than 4 OTUs
sel <- rowSums(otuCompl[, -c(1:2)]) >= 4
otuSel <- otuCompl[sel, ]
envSel <- envCompl[sel, ]


## Removing phosphate
#envSel <- envSel %>% select(-phosphate_mean.tif)


dim(otuSel); dim(envSel)
[1] 983 351
[1] 983 138

Show what samples are left after complete cases and >4 OTUs filters

env_sel_locations<- ggplot(data = envSel,
                      aes(x = X.y,
                          y = Y)) +
  theme_classic() +
  geom_point(aes(colour = bathy),
             size = 2) +
  scale_colour_gradient2(low = "red",
                         mid = "yellow",
                         high = "green") +
  ggtitle("Location of samples in sorted env file")

env_sel_locations

Show what got removed in complete cases filter

envCompl_ccrem<-env_sort%>%filter(!SampID%in%envCompl$SampID)

envCompl_ccrem_locations<- ggplot(data = envCompl_ccrem,
                      aes(x = X.y,
                          y = Y)) +
  theme_classic() +
  geom_point(aes(colour = bathy),
             size = 2) +
  scale_colour_gradient2(low = "red",
                         mid = "yellow",
                         high = "green") +
  ggtitle("Location of removed complete case samples resulting in envSel file")

envCompl_ccrem_locations

Show what got removed in <4 OTUs filter

inv.sel <- rowSums(otuCompl[, -c(1:2)]) < 4
env.invSel <- envCompl[inv.sel, ]

env.invSel_locations<- ggplot(data = env.invSel,
                      aes(x = X.y,
                          y = Y)) +
  theme_classic() +
  geom_point(aes(colour = bathy),
             size = 2) +
  scale_colour_gradient2(low = "red",
                         mid = "yellow",
                         high = "green") +
  ggtitle("Location of removed samples with <4 OTUs resulting in envSel file")

env.invSel_locations

[Optional] Data thinning

# otu_red <- otu1[1:500, ]
# otu <- otu_red
# 
# env_red <- env[1:500, ]
# env <- env_red

Whole dataset for first run

otu<-otuSel[,-c(1)]
env<-envSel[, -c(1,3:23)]

table(is.na(otu))

 FALSE 
344050 
table(is.na(env))

 FALSE 
114028 
#str(otuSel)
#str(envSel)

Splitting in subsets

# ## Class 1
# otu1 <- subset(otuSel, envSel$SplitRev == 1)
# env1 <- envSel %>% filter(SplitRev == 1)
# 
# ## Class 2
# otu2 <- subset(otuSel, envSel$SplitRev == 2)
# env2 <- envSel %>% filter(SplitRev == 2)
# 
# ## Class 4
# otu3 <- subset(otuSel, envSel$SplitRev == 3)
# env3 <- envSel %>% filter(SplitRev == 3)
# 
# ## Class 4
# otu4 <- subset(otuSel, envSel$SplitRev == 4)
# env4 <- envSel %>% filter(SplitRev == 4)
# 
# ## Class 6
# otu6 <- subset(otuSel, envSel$SplitRev == 6)
# env6 <- envSel %>% filter(SplitRev == 6)
# 
# ## Class 7
# otu7 <- subset(otuSel, envSel$SplitRev == 7)
# env7 <- envSel %>% filter(SplitRev == 7)
# 
# ## Class 8
# otu8 <- subset(otuSel, envSel$SplitRev == 8)
# env8 <- envSel %>% filter(SplitRev == 8)

Selecting subset

# ## Selecting 
# otu <- otu2
# env <- env2

format data correctly

#env$coords.x1<-as.numeric(env$coords.x1)
#env$coords.x2<-as.numeric(env$coords.x2)

Abundance weighting

Make function You might have to drop variables that have been imported as character


otu_pa <- decostand(x = otu[, -c(1)],
                    method = "pa")

# y = ax^w         # power transformation formula
dt <- otu[, -c(1:2)]          # species data to transform
x_mn <- min(dt[dt > 0])
x_mx <- max(dt)
rng <- 6           # abundance range
w <- log(rng) / (log(x_mx) - log(x_mn))
a <- x_mn^(-w)
# otu_6 <- a * dt[, -c(1:3)]^w
otu_6 <- a * dt^w
range(otu_6)
[1] 0 6

Ordination methods

DCA PA

Sys.time()
[1] "2022-09-23 17:20:18 CEST"
dca_pa <- decorana(veg = otu_pa)
some species were removed because they were missing in the data
print(dca_pa, head=T)

Call:
decorana(veg = otu_pa) 

Detrended correspondence analysis with 26 segments.
Rescaling of axes with 4 iterations.

                  DCA1   DCA2   DCA3   DCA4
Eigenvalues     0.5071 0.4546 0.2261 0.1962
Decorana values 0.7821 0.4726 0.2896 0.2175
Axis lengths    3.3459 4.4457 4.5141 3.6860

GNMDS PA

Distances - don’t run if loading saved object

## Bray-Curtis
dist_pa <- vegdist(x = otu_pa, method = "bray")

## Geodist
ep <- 0.8     # epsilon
geodist_pa <- isomapdist(dist = dist_pa, epsilon = ep)
Save the result - don’t run if loading saved object
saveRDS(geodist_pa,
        file = (file.path(dataPath,"inputs/LoDens_geodist_pa.rds")))

Ordination

monoMDS - don’t run if loading saved object


Sys.time()
[1] "2022-09-23 17:20:25 CEST"
d <- 2
mds_pa <- list()

for (i in 1:100) {
  mds_pa[[i]]<-monoMDS(geodist_pa,
                    matrix(c(runif(dim(otu_pa)[1]*d)),
                           nrow = dim(otu_pa)[1]),
                    k = d,
                    model = "global",
                    maxit = 2000,
                    smin = 1e-7,
                    sfgrmin = 1e-7)
}
Save the result - don’t run if loading saved object
saveRDS(mds_pa,
        file = file.path(dataPath,"inputs/LoDens_mds_pa.rds"))
Best nmds solution - PA

Make function?

# Loading geodist object
# geodist_nfi <- readRDS(file = "../SplitRev2_geodist_pa_full.rds")
# Loading mds results
# mds <- readRDS(file = "../SplitRev2_mds_pa_full.rds")

## Extracting the stress of each nmds iteration
mds_stress_pa<-unlist(lapply(mds_pa, function(v){v[[22]]})) 

ordered_pa <-order(mds_stress_pa)

## Best, second best, and worst solution
mds_stress_pa[ordered_pa[1]]
[1] 0.2281822
mds_stress_pa[ordered_pa[2]]
[1] 0.2292384
mds_stress_pa[ordered_pa[10]]
[1] 0.2314408
## Scaling of axes to half change units and varimax rotation by postMDS
mds_best_pa<-postMDS(mds_pa[[ordered_pa[1]]],
                  geodist_pa, 
                  pc = TRUE, 
                  halfchange = TRUE, 
                  threshold = ep)     # Is this threshold related to the epsilon above?
mds_best_pa

Call:
monoMDS(dist = geodist_pa, y = matrix(c(runif(dim(otu_pa)[1] *      d)), nrow = dim(otu_pa)[1]), k = d, model = "global", maxit = 2000,      smin = 1e-07, sfgrmin = 1e-07) 

Non-metric Multidimensional Scaling

983 points, dissimilarity ‘bray shortest isomap’, call ‘isomapdist(dist = dist_pa, epsilon = ep)’

Dimensions: 2 
Stress:     0.2281822 
Stress type 1, weak ties
Scores scaled to unit root mean square, rotated to principal components
Stopped after 147 iterations: Stress nearly unchanged (ratio > sratmax)
mds_secbest_pa <- postMDS(mds_pa[[ordered_pa[2]]],
                          geodist_pa, 
                          pc = TRUE, 
                          halfchange = TRUE, 
                          threshold = ep)
mds_secbest_pa

Call:
monoMDS(dist = geodist_pa, y = matrix(c(runif(dim(otu_pa)[1] *      d)), nrow = dim(otu_pa)[1]), k = d, model = "global", maxit = 2000,      smin = 1e-07, sfgrmin = 1e-07) 

Non-metric Multidimensional Scaling

983 points, dissimilarity ‘bray shortest isomap’, call ‘isomapdist(dist = dist_pa, epsilon = ep)’

Dimensions: 2 
Stress:     0.2292384 
Stress type 1, weak ties
Scores scaled to unit root mean square, rotated to principal components
Stopped after 249 iterations: Scale factor of gradient nearly zero (< sfgrmin)
## Procrustes comparisons
procr_pa <- procrustes(mds_best_pa,
                    mds_secbest_pa,
                    permutations=999)
protest(mds_best_pa,
        mds_secbest_pa,
        permutations=999)

Call:
protest(X = mds_best_pa, Y = mds_secbest_pa, permutations = 999) 

Procrustes Sum of Squares (m12 squared):        0.1238 
Correlation in a symmetric Procrustes rotation: 0.9361 
Significance:  0.001 

Permutation: free
Number of permutations: 999
plot(procr_pa)

png(file=file.path(dataPath,"outputs/LoDens_procrustes_pa.png"), width=1000, height=700)
plot(procr_pa)
dev.off()
png 
  2 

Correlation of axis: DCA vs NMDS - PA
# Extracting ordination axis
ax <- 2
axis_pa <- cbind(mds_best_pa$points,
                 scores(dca_pa,
                        display = "sites",
                        origin = TRUE)[, 1:ax])

ggcorr(axis_pa, 
       method=c("everything","kendall"), 
       label = TRUE,
       label_size = 3, 
       label_color = "black",  
       nbreaks = 8,
       label_round = 3,
       low = "red",
       mid = "white",
       high = "green")

NA
NA
Save the figure
ggsave(filename = file.path(dataPath,"outputs/LoDens_correlationPCAvsNMDS_PA.png"),
       device = "png",
       dpi=300 )
Saving 7 x 7 in image
# Switching direction of NMDS1
# mds_best$points[, 1] <- -mds_best$points[, 1]

DCA R6

dca_r6 <- decorana(veg = otu_6)
some species were removed because they were missing in the data
print(dca_r6, head=T)

Call:
decorana(veg = otu_6) 

Detrended correspondence analysis with 26 segments.
Rescaling of axes with 4 iterations.

                  DCA1   DCA2   DCA3   DCA4
Eigenvalues     0.5086 0.4309 0.2327 0.1707
Decorana values 0.8101 0.4269 0.2892 0.2511
Axis lengths    3.9171 4.5861 4.2660 4.4962

GNMDS R6

Distances - don’t run if loading saved object

## Bray-Curtis
dist_r6 <- vegdist(x = otu_6, method = "bray")

## Geodist
ep <- 0.80     # epsilon
geodist_r6 <- isomapdist(dist = dist_r6, epsilon = ep)
Save the result - don’t run if loading saved object
saveRDS(geodist_r6,
        file = (file.path(dataPath,"inputs/LoDens_geodist_r6.rds")))

Ordination

200 reps - don’t run if loading saved object

took 20mins

# monoMDS
d <- 2
mds_r6 <- list()

Sys.time()
[1] "2022-09-23 17:30:43 CEST"
for (i in 1:200) {
  mds_r6[[i]]<-monoMDS(geodist_r6,
                    matrix(c(runif(dim(otu_6)[1]*d)),
                           nrow = dim(otu_6)[1]),
                    k = d,
                    model = "global",
                    maxit = 2000,
                    smin = 1e-7,
                    sfgrmin = 1e-7)
}
Sys.time()
[1] "2022-09-23 17:46:18 CEST"
Save the result - don’t run if loading saved object
saveRDS(mds_r6,
        file = file.path(dataPath,"inputs/LoDens_mds_r6.rds")) 
Best nmds solution r6 200 rep
# Loading geodist object
# geodist_nfi <- readRDS(file = "../SplitRev2_geodist_nfi.rds")

# Loading mds results
# mds <- readRDS(file = "../SplitRev2_mds.rds")

## Extracting the stress of each nmds iteration
mds_stress_r6<-unlist(lapply(mds_r6, function(v){v[[22]]})) 

ordered_r6 <-order(mds_stress_r6)

## Best, second best, and worst solution
mds_stress_r6[ordered_r6[1]]
[1] 0.2129278
mds_stress_r6[ordered_r6[2]]
[1] 0.2139067
mds_stress_r6[ordered_r6[100]]
[1] 0.2172448
## Scaling of axes to half change units and varimax rotation by postMDS
mds_best_r6<-postMDS(mds_r6[[ordered_r6[1]]],
                     geodist_r6, 
                     pc = TRUE, 
                     halfchange = TRUE, 
                     threshold = ep)     # Is this threshold related to the epsilon above?
mds_best_r6

Call:
monoMDS(dist = geodist_r6, y = matrix(c(runif(dim(otu_6)[1] *      d)), nrow = dim(otu_6)[1]), k = d, model = "global", maxit = 2000,      smin = 1e-07, sfgrmin = 1e-07) 

Non-metric Multidimensional Scaling

983 points, dissimilarity ‘bray shortest isomap’, call ‘isomapdist(dist = dist_r6, epsilon = ep)’

Dimensions: 2 
Stress:     0.2129278 
Stress type 1, weak ties
Scores scaled to unit root mean square, rotated to principal components
Stopped after 176 iterations: Stress nearly unchanged (ratio > sratmax)
mds_secbest_r6<-postMDS(mds_r6[[ordered_r6[2]]],
                        geodist_r6, 
                        pc = TRUE, 
                        halfchange = TRUE, 
                        threshold = ep)
mds_secbest_r6

Call:
monoMDS(dist = geodist_r6, y = matrix(c(runif(dim(otu_6)[1] *      d)), nrow = dim(otu_6)[1]), k = d, model = "global", maxit = 2000,      smin = 1e-07, sfgrmin = 1e-07) 

Non-metric Multidimensional Scaling

983 points, dissimilarity ‘bray shortest isomap’, call ‘isomapdist(dist = dist_r6, epsilon = ep)’

Dimensions: 2 
Stress:     0.2139067 
Stress type 1, weak ties
Scores scaled to unit root mean square, rotated to principal components
Stopped after 139 iterations: Scale factor of gradient nearly zero (< sfgrmin)
## Procrustes comparisons
procr_r6 <- procrustes(mds_best_r6,
                       mds_secbest_r6,
                       permutations=999)
protest(mds_best_r6,
        mds_secbest_r6,
        permutations=999)

Call:
protest(X = mds_best_r6, Y = mds_secbest_r6, permutations = 999) 

Procrustes Sum of Squares (m12 squared):        0.1608 
Correlation in a symmetric Procrustes rotation: 0.9161 
Significance:  0.001 

Permutation: free
Number of permutations: 999
plot(procr_r6)

png(file.path(dataPath,"outputs/LoDens_procrustes_r6.png"), width=1000, height=700,) #added 1000
plot(procr_r6)
dev.off()
png 
  2 

#### 1000 reps - don’t run if loading saved object Currently commented out as it gained nothing but took extra time. Can be removed in due course.

Correlation of axis: DCA vs NMDS

retain the 200 rep version

# Extracting ordination axis
ax <- 2
axis_r6 <- cbind(mds_best_r6$points,
                 scores(dca_r6,
                        display = "sites",
                        origin = TRUE)[, 1:ax])

ggcorr(axis_r6, 
       method=c("everything","kendall"), 
       label = TRUE,
       label_size = 3, 
       label_color = "black",  
       nbreaks = 8,
       label_round = 3,
       low = "red",
       mid = "white",
       high = "green")

Save the figure
ggsave(filename = file.path(dataPath,"outputs/LoDens_correlationPCAvsNMDS_r6.png"),
       device = "png",
       dpi=300 )
Saving 7 x 7 in image
# Switching direction of NMDS1
# mds_best$points[, 1] <- -mds_best$points[, 1]

Plotting DCA & GNMDS

## Adding scores to data frame
otu_6$gnmds1 <- mds_best_r6$points[, 1]
otu_6$gnmds2 <- mds_best_r6$points[, 2]
otu_6$dca1 <- scores(dca_r6, display = "sites", origin = TRUE)[, 1]
otu_6$dca2 <- scores(dca_r6, display = "sites", origin = TRUE)[, 2]

p_gnmds_r6 <- ggplot(data = otu_6,
                     aes(x = gnmds1,
                         y = gnmds2)) +
  theme_classic() +
  coord_fixed() +
  ggtitle("GNMDS",
          subtitle = "First run") +
  geom_point(colour = "red") +
  geom_vline(xintercept = 0,
             linetype = 2,
             colour = "lightgray") +
  geom_hline(yintercept = 0,
             linetype = 2,
             colour = "lightgray")

p_dca_r6 <- ggplot(data = otu_6,
                   aes(x = dca1,
                       y = dca2)) +
  theme_classic() +
  coord_fixed() +
  ggtitle("DCA",
          subtitle = "First run") +
  geom_point(colour = "red") +
  geom_vline(xintercept = 0,
             linetype = 2,
             colour = "lightgray") +
  geom_hline(yintercept = 0,
             linetype = 2,
             colour = "lightgray")

p_gnmds_r6 + p_dca_r6

NA
NA
Save the figure
ggsave(file.path(dataPath,"outputs/LoDens_gnmds_dca_r6.png"),
       device = "png", 
       dpi=300)
Saving 7 x 7 in image

Species-environment relationships

Selecting ordination

ord <- mds_best_r6

## Axis scores if selected ord is GNMDS
axis <- ord$points %>% as.data.frame

## Axis scores if selected ord is DCA
# axis <- scores(ord,
#                display = "sites",
#                origin = TRUE)[, 1:ax])

Create additional variables

decided to make MLD-bathy vars

env<-env %>% 
  mutate ("MLDmean_bathy"=MLDmean_Robinson-(bathy*-1),
          "MLDmin_bathy"=MLDmin_Robinson-(bathy*-1),
          "MLDmax_bathy"=MLDmax_Robinson-(bathy*-1))

env$MLDmean_bathy<-cut(env$MLDmean_bathy, 
      breaks=c(-2560, -20,20,130),#checked range of values first (min -2554, max 123)
      labels=c('belowMLD','onPycno','inMixLayer'))
env$MLDmin_bathy<-cut(env$MLDmin_bathy, 
      breaks=c(-2560, -20,20,130),#checked range of values first (min -2554, max 123)
      labels=c('belowMLD','onPycno','inMixLayer'))
env$MLDmax_bathy<-cut(env$MLDmax_bathy, 
      breaks=c(-2560, -20,20,130),#checked range of values first (min -2554, max 123)
      labels=c('belowMLD','onPycno','inMixLayer'))

env$swDensRob_avs<-swRho(salinity=env$Smean_Robinson,
                         temperature=env$Tmean_Robinson,
                         pressure=(env$bathy*-1),
                         eos="unesco")

Correlation ordination axes and environmental variables

Removing non-env vars

env_cont<-env%>% select(-c(landscape,sedclass,gmorph, MLDmean_bathy, MLDmax_bathy, MLDmin_bathy, #categorical
                           optional, #not a var
                           MLDmax_Robinson, MLDmean_Robinson, MLDmin_Robinson,  #replaced by new vars
                           MLDsd_Robinson #not meaningful
                           ))
env_cont<-env_cont%>% mutate_if(is.integer,as.numeric)
env_corr <- env_cont %>% select(-c(SampID))

# env_corr$coords.x1<-as.numeric(env_corr$coords.x1)
# env_corr$coords.x2<-as.numeric(env_corr$coords.x2)

env_corr[(!is.numeric(env_corr)),]

Correlations

# Vector to hold correlations
cor_ax1 <- NULL
cor_ax2 <- NULL
pv_ax1 <- NULL
pv_ax2 <- NULL

# NMDS1
for( i in seq(length(env_corr))) {
  ct.i <- cor.test(axis$MDS1,
                   env_corr[, i],
                   method = "kendall")
  cor_ax1[i] <- ct.i$estimate
  pv_ax1[i] <- ct.i$p.value
}

# NMDS2
for( i in seq(length(env_corr))) {
  ct.i <- cor.test(axis$MDS2,
                   env_corr[, i],
                   method = "kendall")
  cor_ax2[i] <- ct.i$estimate
  pv_ax2[i] <- ct.i$p.value
}

cor_tab <- data.frame(env = names(env_corr),
                      ord_ax1 = cor_ax1,
                      pval_ax1 = pv_ax1,
                      ord_ax2 = cor_ax2,
                      pval_ax2 = pv_ax2)

cor_tab

write.csv(x = cor_tab,
          file = file.path(dataPath,"inputs/LoDens_cor-table_r6_200rep_MLD-bathy.csv"),
          row.names = FALSE)

Dot chart to check for gaps in correlation

cor_a1_sort<-cor_tab%>%
  mutate(abs_ord_ax1=abs(ord_ax1),
         abs_ord_ax2=abs(ord_ax2)) %>%
  arrange(desc(abs_ord_ax1))

cor_a2_sort<-cor_tab%>%
  mutate(abs_ord_ax1=abs(ord_ax1),
         abs_ord_ax2=abs(ord_ax2)) %>%
  arrange(desc(abs_ord_ax2))

dotchart(cor_a1_sort$abs_ord_ax1, main="Absolute (+/-) correlations between envVars and gnmds axis 1")


cor_cut<-0.37 #decide

cor_sel<-subset(cor_a1_sort,abs_ord_ax1>cor_cut)
cor_sel

Sel env var (top corr)

env_os <- env[, cor_sel$env]
env_os
str(env_os)
'data.frame':   983 obs. of  11 variables:
 $ BO22_icecoverrange_ss  : num  0 0 0 0 0 0 0 0 0 0 ...
 $ BO22_icethickrange_ss  : num  0 0 0 0 0 0 0 0 0 0 ...
 $ BO22_icethickltmax_ss  : num  0 0 0 0 0 0 0 0 0 0 ...
 $ BO22_icecoverltmax_ss  : num  0 0 0 0 0 0 0 0 0 0 ...
 $ BO22_icethickmean_ss   : num  0 0 0 0 0 0 0 0 0 0 ...
 $ BO22_icecovermean_ss   : num  0 0 0 0 0 0 0 0 0 0 ...
 $ Tmean_Robinson         : num  7.87 7.85 6.92 4.27 7.36 ...
 $ Tmax_Robinson          : num  8.03 8.03 7.08 4.62 7.44 ...
 $ Tmin_Robinson          : num  7.72 7.74 6.6 3.92 7.23 ...
 $ temp_mean              : num  6.22 5.93 4.33 1.78 5.34 ...
 $ BO22_dissoxltmin_bdmean: num  262 263 263 264 263 ...

Ordisurfs top corr

ordsrfs <- list(length = ncol(env_os))

for (i in seq(ncol(env_os))) {
  os.i <- gg_ordisurf(ord = ord,
                      env.var = env_os[, i],
                      pt.size = 1,
                      # binwidth = 0.05,
                      var.label = names(env_os)[i],
                      gen.text.size = 10,
                      title.text.size = 15,
                      leg.text.size = 10)
  
  ordsrfs[[i]] <- os.i$plot
}


ordsrfs_plt <- ggarrange(plotlist = ordsrfs,
                         nrow = 4,
                         ncol = 3)

ordsrfs_plt

Save some outputs
ggexport(ordsrfs_plt,
          filename = file.path(dataPath,"outputs/LoDens_ordisurfs_top_corr.png"),
          width = 2000,
          height = 2500)

Sel env var (manual)

env_os_m <- env[,c("Tmean_Robinson", #top corr
                  "salt_max", #top corr
                  "Smax_Robinson", #comparison to top corr
                  "swDensRob_avs", #top corr
                  "BO22_icecoverltmax_ss",#top corr ax2
                  "BO22_icecovermean_ss",#top corr
                  "BO22_dissoxmean_bdmean",#top corr
                  #"BO22_carbonphytoltmin_bdmean",#top corr - no clear gradient in ordisurf
                  "BO22_ppltmin_ss", #top corr
                  "X.y", #comparison to Y
                  "Y", #top corr
                  "spd_std", #top corr ax2 (blended model)
                  "CSpdsd_Robinson", #comparison to top corr ax2 (blended model)
                  "mud", #highest sed var ax1 + corr
                  "gravel",#highest sed var ax1 - corr
                  "BO22_silicateltmax_bdmean", #just under top corr ax1
                  "bathy" #intuitive for comparisons
                )]
env_os_m
str(env_os_m)
'data.frame':   983 obs. of  16 variables:
 $ Tmean_Robinson           : num  7.87 7.85 6.92 4.27 7.36 ...
 $ salt_max                 : num  35.2 35.2 35.2 35.1 35.2 ...
 $ Smax_Robinson            : num  35.4 35.4 35.4 35.2 35.4 ...
 $ swDensRob_avs            : num  1029 1029 1030 1031 1030 ...
 $ BO22_icecoverltmax_ss    : num  0 0 0 0 0 0 0 0 0 0 ...
 $ BO22_icecovermean_ss     : num  0 0 0 0 0 0 0 0 0 0 ...
 $ BO22_dissoxmean_bdmean   : num  285 286 286 286 286 ...
 $ BO22_ppltmin_ss          : num  0.000174 0.000171 0.00016 0.000153 0.000182 ...
 $ X.y                      : num  76461 67461 57261 50461 61261 ...
 $ Y                        : num  7200734 7206934 7209334 7203134 7191534 ...
 $ spd_std                  : num  0.074 0.096 0.114 0.098 0.114 ...
 $ CSpdsd_Robinson          : num  0.01692 0.00818 0.0035 0.00648 0.00254 ...
 $ mud                      : num  26 18 8 8 1 26 18 1 7 1 ...
 $ gravel                   : num  36 55 21 21 33 36 55 35 65 35 ...
 $ BO22_silicateltmax_bdmean: num  6.29 6.27 7.41 9.03 6.91 ...
 $ bathy                    : num  -325 -334 -481 -604 -435 ...

Ordisurfs top corr

ordsrfs_m <- list(length = ncol(env_os_m))

for (i in seq(ncol(env_os_m))) {
  os.i_m <- gg_ordisurf(ord = ord,
                      env.var = env_os_m[, i],
                      pt.size = 1,
                      # binwidth = 0.05,
                      var.label = names(env_os_m)[i],
                      gen.text.size = 10,
                      title.text.size = 15,
                      leg.text.size = 10)
  
  ordsrfs_m[[i]] <- os.i_m$plot
}


ordsrfs_plt_m <- ggarrange(plotlist = ordsrfs_m,
                         nrow = 4,
                         ncol = 4)

ordsrfs_plt_m

Save some outputs
ggexport(ordsrfs_plt_m,
          filename = file.path(dataPath,"outputs/LoDens_ordisurfs_man_sel.png"),
          width = 2000,
          height = 2000)

Envfit

## Select if any var should be excluded from envfit (makes less busy to read)
env_os_m_envfit<-env_os_m [,c("Tmean_Robinson", #top corr
                  "salt_max", #top corr
                  "Smax_Robinson", #comparison to top corr
                  "swDensRob_avs", #top corr
                  "BO22_icecoverltmax_ss",#top corr ax2
                  "BO22_icecovermean_ss",#top corr
                  "BO22_dissoxmean_bdmean",#top corr
                  #"BO22_carbonphytoltmin_bdmean",#top corr - no clear gradient in ordisurf
                  "BO22_ppltmin_ss", #top corr
                  "X.y", #comparison to Y
                  "Y", #top corr
                 # "spd_std", #top corr ax2 (blended model)
                  "CSpdsd_Robinson", #comparison to top corr ax2 (blended model)
                  "mud", #highest sed var ax1 + corr
                  "gravel",#highest sed var ax1 - corr
                  "BO22_silicateltmax_bdmean", #just under top corr ax1
                  "bathy" #intuitive for comparisons
                
  
)]

colnames(env_os_m_envfit)<-c("T", #top corr
                  "sMx____", #top corr
                  "SmaxR", #comparison to top corr
                  "swDensR", #top corr
                  "icecovmax",#top corr ax2
                  "icecovav",#top corr
                  "dissoxav",#top corr
                  #"BO22_carbonphytoltmin_bdmean",#top corr - no clear gradient in ordisurf
                  "ppltmin", #top corr
                  "X", #comparison to Y
                  "Y", #top corr
                 # "spd_std", #top corr ax2 (blended model)
                  "CSsd", #comparison to top corr ax2 (blended model)
                  "mud", #highest sed var ax1 + corr
                  "gravel",#highest sed var ax1 - corr
                  "SiLtmax", #just under top corr ax1
                  "bathy" #intuitive for comparisons
                 )

## Envfot plot
gg_envfit(ord = ord,
          env = env_os_m_envfit,
          pt.size = 1)


## Envfit analysis

ef <- envfit(ord = ord,
             env = env_os_m_envfit,
             # na.rm = TRUE
             )
efDF <- as.data.frame(scores(ef,
                             display = "vectors"))
Save the plot
ggsave(filename = file.path(dataPath,"outputs/LoDens_EnvFit_man_sel_cln.png"),
       device = "png",
       dpi=300 )
Saving 7 x 7 in image

Categorical envVar visualised on the mdsplots

use dataset with categorical var included

env_vis<-env
env_vis$gnmds1 <- otu_6$gnmds1
env_vis$gnmds2 <- otu_6$gnmds2

gnmds w mld mean - bathy

p_mld <- ggplot(data = env_vis,
                     aes(x = gnmds1,
                         y = gnmds2)) +
  theme_classic() +
  coord_fixed() +
  ggtitle("GNMDS coloured by proximity to mixed layer depth",
          subtitle = "First run") +
  geom_point(aes(colour = MLDmean_bathy)) +
  geom_vline(xintercept = 0,
             linetype = 2,
             colour = "lightgray") +
  geom_hline(yintercept = 0,
             linetype = 2,
             colour = "lightgray")

p_mld

Save the plot
ggsave(filename = file.path(dataPath,"outputs/LoDens_gnmds_MLDmeanBathy.png"),
       device = "png",
       dpi=300 )
Saving 7 x 7 in image

gnmds w sedclass

deal with sedclass codes
env_vis<- env_vis %>% 
  mutate(
    sedclassName = case_when(
            sedclass == "1" ~ "SedCoverR",
            sedclass == "5" ~ "Rock",
            sedclass == "20" ~ "Mud",
            sedclass == "21" ~ "MwBlock",
            sedclass == "40" ~ "sMud",
            sedclass == "80" ~ "mSand",
            sedclass == "100" ~ "Sand",
            sedclass == "110" ~ "gMud",
            sedclass == "115" ~ "gsMud",
            sedclass == "120" ~ "gmSand",
            sedclass == "130" ~ "gSand",
            sedclass == "150" ~ "MSG",
            sedclass == "160" ~ "sGravel",
            sedclass == "170" ~ "Gravel",
            sedclass == "175" ~ "GravBlock",
            sedclass == "185" ~ "SGBmix",
            sedclass == "205" ~ "S/MwB",
            sedclass == "206" ~ "S/MwG/B",
            sedclass == "215" ~ "SGBalt",
            sedclass == "300" ~ "HardSed",
            sedclass == "500" ~ "Biogenic"
                        
    )
  )
colour palette to cope with up to 25 categorical colours
c25 <- c(
  "dodgerblue2", "#E31A1C", # red
  "green4",
  "#6A3D9A", # purple
  "#FF7F00", # orange
  "black", "gold1",
  "skyblue2", "#FB9A99", # lt pink
  "palegreen2",
  "#CAB2D6", # lt purple
  "#FDBF6F", # lt orange
  "gray70", "khaki2",
  "maroon", "orchid1", "deeppink1", "blue1", "steelblue4",
  "darkturquoise", "green1", "yellow4", "yellow3",
  "darkorange4", "brown"
)

p_sed <- ggplot(data = env_vis,
                     aes(x = gnmds1,
                         y = gnmds2)) +
  theme_classic() +
  coord_fixed() +
  ggtitle("GNMDS coloured by sediment class",
          subtitle = "First run") +
  geom_point(aes(colour = factor(sedclassName))) +
   scale_colour_manual(values=c25)+
  geom_vline(xintercept = 0,
             linetype = 2,
             colour = "lightgray") +
  geom_hline(yintercept = 0,
             linetype = 2,
             colour = "lightgray")+
  guides(colour=guide_legend(ncol=2))

p_sed

Save the plot
ggsave(filename = file.path(dataPath,"outputs/LoDens_gnmds_sedclass.png"),
       device = "png",
       dpi=300 )
Saving 7 x 7 in image

gnmds w lanscape

deal with sedclass codes
env_vis<- env_vis %>% 
  mutate(
    landscapeName = case_when(
            landscape == "1" ~ "Strandflat",
            landscape == "21" ~ "ContSlope",
            landscape == "22" ~ "Canyon",
            landscape == "31" ~ "Valley",
            landscape == "32" ~ "Fjord",
            landscape == "41" ~ "DeepPlain",
            landscape == "42" ~ "SlopePlain",
            landscape == "43" ~ "ShelfPlain",
            landscape == "431" ~ "shallowValley"
    )
  )

p_land <- ggplot(data = env_vis,
                     aes(x = gnmds1,
                         y = gnmds2)) +
  theme_classic() +
  coord_fixed() +
  ggtitle("GNMDS coloured by landscape class",
          subtitle = "First run") +
  geom_point(aes(colour = factor(landscapeName))) +
  geom_vline(xintercept = 0,
             linetype = 2,
             colour = "lightgray") +
  geom_hline(yintercept = 0,
             linetype = 2,
             colour = "lightgray")+
  guides(colour=guide_legend(ncol=2))

p_land

Save the plot
ggsave(filename = file.path(dataPath,"outputs/LoDens_gnmds_landscape.png"),
       device = "png",
       dpi=300 )
Saving 7 x 7 in image

Gnmds w gmorph


p_gmo <- ggplot(data = env_vis,
                     aes(x = gnmds1,
                         y = gnmds2)) +
  theme_classic() +
  coord_fixed() +
  ggtitle("GNMDS coloured by landscape class",
          subtitle = "First run") +
  geom_point(aes(colour = factor(gmorph))) +
  geom_vline(xintercept = 0,
             linetype = 2,
             colour = "lightgray") +
  geom_hline(yintercept = 0,
             linetype = 2,
             colour = "lightgray")+
  guides(colour=guide_legend(ncol=2))

p_gmo

Save the plot
ggsave(filename = file.path(dataPath,"outputs/LoDens_gnmds_gmorph.png"),
       device = "png",
       dpi=300 )
Saving 7 x 7 in image
cat_var_plots<-p_mld+p_sed+p_land+p_gmo
Save the plot
ggexport(cat_var_plots,
          filename = file.path(dataPath,"outputs/LoDens_gnmds_catvar.png"),
          width = 1000,
          height = 800)

Sample identification in the mdsplot

p_gmo <- ggplot(data = env_vis,
                     aes(x = gnmds1,
                         y = gnmds2)) +
  theme_classic() +
  coord_fixed() +
  ggtitle("GNMDS coloured by landscape class",
          subtitle = "First run") +
  geom_point(aes(colour = factor(SampID))) +
  geom_vline(xintercept = 0,
             linetype = 2,
             colour = "lightgray") +
  geom_hline(yintercept = 0,
             linetype = 2,
             colour = "lightgray")+
  guides(colour="none", size="none")

ggplotly(p_gmo)
LS0tDQp0aXRsZTogIk1BUkVBTk8gLSBOaU4gLSBsb0RlbnMgPDEwMzEuOCINCmF1dGhvcnM6IFRoaWpzIHZhbiBTb24sIFJ1bmUgSGFsdm9yc2VuLCBSZWJlY2NhIFJvc3MsIEdlbm92ZXZhIEdvbnphbGVzLU1pcmVsaXMsIE1hcmdhcmV0IERvbGFuDQpkYXRlOiAiTGFzdCBSZW5kZXJlZCBvbiBgciBmb3JtYXQoU3lzLnRpbWUoKSwgJyVkICVCLCAlWScpYCINCm91dHB1dDogDQogIGh0bWxfbm90ZWJvb2s6IA0KICAgIHRvYzogeWVzDQogICAgdG9jX2RlcHRoOiAyDQogICAgdG9jX2Zsb2F0OiB5ZXMNCiAgICBmaWdfd2lkdGg6IDcNCiAgICBmaWdfaGVpZ2h0OiA3DQphbHdheXNfYWxsb3dfaHRtbDogdHJ1ZSANCi0tLQ0KDQojIyBQYWNrYWdlcw0KYGBge3J9DQpsaWJyYXJ5KGRwbHlyKQ0KbGlicmFyeShnZ2NvcnJwbG90KQ0KbGlicmFyeShnZ3B1YnIpDQpsaWJyYXJ5KHBhdGNod29yaykNCiMgbGlicmFyeSgiaGVyZSIpDQojIGxpYnJhcnkoImJvb2tkb3duIikNCiMgbGlicmFyeSgiZG93bmxvYWR0aGlzIikNCmxpYnJhcnkodmVnYW4pDQpsaWJyYXJ5KHBseXIpDQpsaWJyYXJ5KGUxMDcxKQ0KbGlicmFyeSh0aWR5dmVyc2UpDQpsaWJyYXJ5KHZpcmlkaXMpDQpsaWJyYXJ5KEdHYWxseSkNCmxpYnJhcnkoZ2dyZXBlbCkNCmxpYnJhcnkoZ2dvcmRpcGxvdHMpDQpsaWJyYXJ5KHJlYWRyKQ0KbGlicmFyeShSQ29sb3JCcmV3ZXIpDQpsaWJyYXJ5KG9jZSkgDQpsaWJyYXJ5KHBsb3RseSkNCg0KDQpzb3VyY2UoImdnX29yZGlzdXJmX3ZpcmlkaXMuUiIpDQpgYGANCg0KDQoNCiMjIERhdGENCmBgYHtyfQ0KDQplbnZfc29ydCA8LSByZWFkLmNzdihmaWxlLnBhdGgoZGF0YVBhdGgsImlucHV0cy9MRGVuc1NobHdfZW52X3NvcnRfMjAyMi0wOS0yMy5jc3YiKSkgJT4lIGFzLmRhdGEuZnJhbWUNCm90dV9zb3J0IDwtIHJlYWQuY3N2KGZpbGUucGF0aChkYXRhUGF0aCwiaW5wdXRzL0xEZW5zU2hsd19vdHVfc29ydF8yMDIyLTA5LTIzLmNzdiIpKSAlPiUgYXMuZGF0YS5mcmFtZQ0KDQoNCnRhYmxlKGVudl9zb3J0JFNhbXBJRDIgPT0gb3R1X3NvcnQkU2FtcElEKQ0KDQoNCiNtdXN0IGJlIHNhbWUgbGVuZ3RoIGFuZCBlcXVhbCB0byB1bmlxdWUgc2FtcGxlIG51bWJlciBsZW5ndGggLSB0aGUgb3JkZXJlZCBsaXN0cyBhcmUgYXNzdW1lZCB0byBiZSBkaXJlY3RseSByZWxhdGFibGUgb24gYSByb3cgYnkgcm93IGJhc2lzDQoNCmRpbShlbnZfc29ydCkNCmRpbShvdHVfc29ydCkNCmBgYA0KDQojIyMjIHJlbW92ZSBhbnkgdmFyaWFibGVzIHdpdGggbm90IGVub3VnaCBjb3ZlcmFnZQ0KYGBge3J9DQoNCmVudl9zb3J0IDwtIGVudl9zb3J0ICU+JSBzZWxlY3QoLWMoIkJPMjJfbGlnaHRib3RtZWFuX2JkbWVhbiIsDQogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIkJPMjJfbGlnaHRib3RsdG1heF9iZG1lYW4iLA0KICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICJCTzIyX2xpZ2h0Ym90bHRtaW5fYmRtZWFuIiwNCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiQk8yMl9saWdodGJvdHJhbmdlX2JkbWVhbiIpKQ0KYGBgDQoNCg0KIyMjIyBQbG90IG1hcCB0byBjaGVjayBsb2NhdGlvbiBvZiBzYW1wbGVzDQpgYGB7cn0NCmVudl9zb3J0X2xvY2F0aW9uczwtIGdncGxvdChkYXRhID0gZW52X3NvcnQsDQogICAgICAgICAgICAgICAgICAgICAgYWVzKHggPSBYLnksDQogICAgICAgICAgICAgICAgICAgICAgICAgIHkgPSBZKSkgKw0KICB0aGVtZV9jbGFzc2ljKCkgKw0KICBnZW9tX3BvaW50KGFlcyhjb2xvdXIgPSBiYXRoeSksDQogICAgICAgICAgICAgc2l6ZSA9IDIpICsNCiAgc2NhbGVfY29sb3VyX2dyYWRpZW50Mihsb3cgPSAicmVkIiwNCiAgICAgICAgICAgICAgICAgICAgICAgICBtaWQgPSAieWVsbG93IiwNCiAgICAgICAgICAgICAgICAgICAgICAgICBoaWdoID0gImdyZWVuIikgKw0KICBnZ3RpdGxlKCJMb2NhdGlvbiBvZiBzYW1wbGVzIGluIHNvcnRlZCBlbnYgZmlsZSIpDQoNCmVudl9zb3J0X2xvY2F0aW9ucw0KYGBgDQoNCg0KIyMjIERhdGEgY2xlYW5pbmcNCmBgYHtyfQ0KIyMgUmVtb3ZpbmcgTkFzDQpvdHVDb21wbCA8LSBvdHVfc29ydFtjb21wbGV0ZS5jYXNlcyhlbnZfc29ydFssIC1jKDE6MjMpXSksIF0NCmVudkNvbXBsIDwtIGVudl9zb3J0W2NvbXBsZXRlLmNhc2VzKGVudl9zb3J0WywgLWMoMToyMyldKSwgXQ0KDQojIyBSZW1vdmluZyBvYnNlcnZhdGlvbnMgd2l0aCBsZXNzIHRoYW4gNCBPVFVzDQpzZWwgPC0gcm93U3VtcyhvdHVDb21wbFssIC1jKDE6MildKSA+PSA0DQpvdHVTZWwgPC0gb3R1Q29tcGxbc2VsLCBdDQplbnZTZWwgPC0gZW52Q29tcGxbc2VsLCBdDQoNCg0KIyMgUmVtb3ZpbmcgcGhvc3BoYXRlDQojZW52U2VsIDwtIGVudlNlbCAlPiUgc2VsZWN0KC1waG9zcGhhdGVfbWVhbi50aWYpDQoNCg0KZGltKG90dVNlbCk7IGRpbShlbnZTZWwpDQpgYGANCiMjIyMgU2hvdyB3aGF0IHNhbXBsZXMgYXJlIGxlZnQgYWZ0ZXIgY29tcGxldGUgY2FzZXMgYW5kID40IE9UVXMgZmlsdGVycw0KYGBge3J9DQplbnZfc2VsX2xvY2F0aW9uczwtIGdncGxvdChkYXRhID0gZW52U2VsLA0KICAgICAgICAgICAgICAgICAgICAgIGFlcyh4ID0gWC55LA0KICAgICAgICAgICAgICAgICAgICAgICAgICB5ID0gWSkpICsNCiAgdGhlbWVfY2xhc3NpYygpICsNCiAgZ2VvbV9wb2ludChhZXMoY29sb3VyID0gYmF0aHkpLA0KICAgICAgICAgICAgIHNpemUgPSAyKSArDQogIHNjYWxlX2NvbG91cl9ncmFkaWVudDIobG93ID0gInJlZCIsDQogICAgICAgICAgICAgICAgICAgICAgICAgbWlkID0gInllbGxvdyIsDQogICAgICAgICAgICAgICAgICAgICAgICAgaGlnaCA9ICJncmVlbiIpICsNCiAgZ2d0aXRsZSgiTG9jYXRpb24gb2Ygc2FtcGxlcyBpbiBzb3J0ZWQgZW52IGZpbGUiKQ0KDQplbnZfc2VsX2xvY2F0aW9ucw0KYGBgDQoNCiMjIyMgU2hvdyB3aGF0IGdvdCByZW1vdmVkIGluIGNvbXBsZXRlIGNhc2VzIGZpbHRlcg0KYGBge3J9DQplbnZDb21wbF9jY3JlbTwtZW52X3NvcnQlPiVmaWx0ZXIoIVNhbXBJRCVpbiVlbnZDb21wbCRTYW1wSUQpDQoNCmVudkNvbXBsX2NjcmVtX2xvY2F0aW9uczwtIGdncGxvdChkYXRhID0gZW52Q29tcGxfY2NyZW0sDQogICAgICAgICAgICAgICAgICAgICAgYWVzKHggPSBYLnksDQogICAgICAgICAgICAgICAgICAgICAgICAgIHkgPSBZKSkgKw0KICB0aGVtZV9jbGFzc2ljKCkgKw0KICBnZW9tX3BvaW50KGFlcyhjb2xvdXIgPSBiYXRoeSksDQogICAgICAgICAgICAgc2l6ZSA9IDIpICsNCiAgc2NhbGVfY29sb3VyX2dyYWRpZW50Mihsb3cgPSAicmVkIiwNCiAgICAgICAgICAgICAgICAgICAgICAgICBtaWQgPSAieWVsbG93IiwNCiAgICAgICAgICAgICAgICAgICAgICAgICBoaWdoID0gImdyZWVuIikgKw0KICBnZ3RpdGxlKCJMb2NhdGlvbiBvZiByZW1vdmVkIGNvbXBsZXRlIGNhc2Ugc2FtcGxlcyByZXN1bHRpbmcgaW4gZW52U2VsIGZpbGUiKQ0KDQplbnZDb21wbF9jY3JlbV9sb2NhdGlvbnMNCmBgYA0KIyMjIyBTaG93IHdoYXQgZ290IHJlbW92ZWQgaW4gPDQgT1RVcyBmaWx0ZXINCg0KDQpgYGB7cn0NCmludi5zZWwgPC0gcm93U3VtcyhvdHVDb21wbFssIC1jKDE6MildKSA8IDQNCmVudi5pbnZTZWwgPC0gZW52Q29tcGxbaW52LnNlbCwgXQ0KDQplbnYuaW52U2VsX2xvY2F0aW9uczwtIGdncGxvdChkYXRhID0gZW52LmludlNlbCwNCiAgICAgICAgICAgICAgICAgICAgICBhZXMoeCA9IFgueSwNCiAgICAgICAgICAgICAgICAgICAgICAgICAgeSA9IFkpKSArDQogIHRoZW1lX2NsYXNzaWMoKSArDQogIGdlb21fcG9pbnQoYWVzKGNvbG91ciA9IGJhdGh5KSwNCiAgICAgICAgICAgICBzaXplID0gMikgKw0KICBzY2FsZV9jb2xvdXJfZ3JhZGllbnQyKGxvdyA9ICJyZWQiLA0KICAgICAgICAgICAgICAgICAgICAgICAgIG1pZCA9ICJ5ZWxsb3ciLA0KICAgICAgICAgICAgICAgICAgICAgICAgIGhpZ2ggPSAiZ3JlZW4iKSArDQogIGdndGl0bGUoIkxvY2F0aW9uIG9mIHJlbW92ZWQgc2FtcGxlcyB3aXRoIDw0IE9UVXMgcmVzdWx0aW5nIGluIGVudlNlbCBmaWxlIikNCg0KZW52LmludlNlbF9sb2NhdGlvbnMNCg0KYGBgDQoNCiMjIyBbT3B0aW9uYWxdIERhdGEgdGhpbm5pbmcNCmBgYHtyLCBldmFsPUZBTFNFfQ0KIyBvdHVfcmVkIDwtIG90dTFbMTo1MDAsIF0NCiMgb3R1IDwtIG90dV9yZWQNCiMgDQojIGVudl9yZWQgPC0gZW52WzE6NTAwLCBdDQojIGVudiA8LSBlbnZfcmVkDQpgYGANCg0KIyMjIFdob2xlIGRhdGFzZXQgZm9yIGZpcnN0IHJ1bg0KYGBge3J9DQpvdHU8LW90dVNlbFssLWMoMSldDQplbnY8LWVudlNlbFssIC1jKDEsMzoyMyldDQoNCnRhYmxlKGlzLm5hKG90dSkpDQoNCnRhYmxlKGlzLm5hKGVudikpDQoNCiNzdHIob3R1U2VsKQ0KI3N0cihlbnZTZWwpDQpgYGANCg0KDQojIyMgU3BsaXR0aW5nIGluIHN1YnNldHMNCmBgYHtyfQ0KIyAjIyBDbGFzcyAxDQojIG90dTEgPC0gc3Vic2V0KG90dVNlbCwgZW52U2VsJFNwbGl0UmV2ID09IDEpDQojIGVudjEgPC0gZW52U2VsICU+JSBmaWx0ZXIoU3BsaXRSZXYgPT0gMSkNCiMgDQojICMjIENsYXNzIDINCiMgb3R1MiA8LSBzdWJzZXQob3R1U2VsLCBlbnZTZWwkU3BsaXRSZXYgPT0gMikNCiMgZW52MiA8LSBlbnZTZWwgJT4lIGZpbHRlcihTcGxpdFJldiA9PSAyKQ0KIyANCiMgIyMgQ2xhc3MgNA0KIyBvdHUzIDwtIHN1YnNldChvdHVTZWwsIGVudlNlbCRTcGxpdFJldiA9PSAzKQ0KIyBlbnYzIDwtIGVudlNlbCAlPiUgZmlsdGVyKFNwbGl0UmV2ID09IDMpDQojIA0KIyAjIyBDbGFzcyA0DQojIG90dTQgPC0gc3Vic2V0KG90dVNlbCwgZW52U2VsJFNwbGl0UmV2ID09IDQpDQojIGVudjQgPC0gZW52U2VsICU+JSBmaWx0ZXIoU3BsaXRSZXYgPT0gNCkNCiMgDQojICMjIENsYXNzIDYNCiMgb3R1NiA8LSBzdWJzZXQob3R1U2VsLCBlbnZTZWwkU3BsaXRSZXYgPT0gNikNCiMgZW52NiA8LSBlbnZTZWwgJT4lIGZpbHRlcihTcGxpdFJldiA9PSA2KQ0KIyANCiMgIyMgQ2xhc3MgNw0KIyBvdHU3IDwtIHN1YnNldChvdHVTZWwsIGVudlNlbCRTcGxpdFJldiA9PSA3KQ0KIyBlbnY3IDwtIGVudlNlbCAlPiUgZmlsdGVyKFNwbGl0UmV2ID09IDcpDQojIA0KIyAjIyBDbGFzcyA4DQojIG90dTggPC0gc3Vic2V0KG90dVNlbCwgZW52U2VsJFNwbGl0UmV2ID09IDgpDQojIGVudjggPC0gZW52U2VsICU+JSBmaWx0ZXIoU3BsaXRSZXYgPT0gOCkNCmBgYA0KDQojIyBTZWxlY3Rpbmcgc3Vic2V0DQpgYGB7cn0NCiMgIyMgU2VsZWN0aW5nIA0KIyBvdHUgPC0gb3R1Mg0KIyBlbnYgPC0gZW52Mg0KDQpgYGANCg0KIyMgZm9ybWF0IGRhdGEgY29ycmVjdGx5DQpgYGB7cn0NCiNlbnYkY29vcmRzLngxPC1hcy5udW1lcmljKGVudiRjb29yZHMueDEpDQojZW52JGNvb3Jkcy54MjwtYXMubnVtZXJpYyhlbnYkY29vcmRzLngyKQ0KDQoNCmBgYA0KDQoNCiMjIEFidW5kYW5jZSB3ZWlnaHRpbmcNCk1ha2UgZnVuY3Rpb24NCllvdSBtaWdodCBoYXZlIHRvIGRyb3AgdmFyaWFibGVzIHRoYXQgaGF2ZSBiZWVuIGltcG9ydGVkIGFzIGNoYXJhY3Rlcg0KYGBge3J9DQoNCm90dV9wYSA8LSBkZWNvc3RhbmQoeCA9IG90dVssIC1jKDEpXSwNCiAgICAgICAgICAgICAgICAgICAgbWV0aG9kID0gInBhIikNCg0KIyB5ID0gYXhedyAgICAgICAgICMgcG93ZXIgdHJhbnNmb3JtYXRpb24gZm9ybXVsYQ0KZHQgPC0gb3R1WywgLWMoMToyKV0gICAgICAgICAgIyBzcGVjaWVzIGRhdGEgdG8gdHJhbnNmb3JtDQp4X21uIDwtIG1pbihkdFtkdCA+IDBdKQ0KeF9teCA8LSBtYXgoZHQpDQpybmcgPC0gNiAgICAgICAgICAgIyBhYnVuZGFuY2UgcmFuZ2UNCncgPC0gbG9nKHJuZykgLyAobG9nKHhfbXgpIC0gbG9nKHhfbW4pKQ0KYSA8LSB4X21uXigtdykNCiMgb3R1XzYgPC0gYSAqIGR0WywgLWMoMTozKV1edw0Kb3R1XzYgPC0gYSAqIGR0XncNCnJhbmdlKG90dV82KQ0KYGBgDQoNCiMjIE9yZGluYXRpb24gbWV0aG9kcw0KIyMjIERDQSBQQQ0KYGBge3J9DQpTeXMudGltZSgpDQpkY2FfcGEgPC0gZGVjb3JhbmEodmVnID0gb3R1X3BhKQ0KDQpwcmludChkY2FfcGEsIGhlYWQ9VCkNCmBgYA0KDQoNCiMjIyBHTk1EUyBQQQ0KDQojIyMjIERpc3RhbmNlcyAtIGRvbid0IHJ1biBpZiBsb2FkaW5nIHNhdmVkIG9iamVjdA0KDQpgYGB7cn0NCiMjIEJyYXktQ3VydGlzDQpkaXN0X3BhIDwtIHZlZ2Rpc3QoeCA9IG90dV9wYSwgbWV0aG9kID0gImJyYXkiKQ0KDQojIyBHZW9kaXN0DQplcCA8LSAwLjggICAgICMgZXBzaWxvbg0KZ2VvZGlzdF9wYSA8LSBpc29tYXBkaXN0KGRpc3QgPSBkaXN0X3BhLCBlcHNpbG9uID0gZXApDQpgYGANCg0KDQojIyMjIyBTYXZlIHRoZSByZXN1bHQgLSBkb24ndCBydW4gaWYgbG9hZGluZyBzYXZlZCBvYmplY3QNCmBgYHtyfQ0Kc2F2ZVJEUyhnZW9kaXN0X3BhLA0KICAgICAgICBmaWxlID0gKGZpbGUucGF0aChkYXRhUGF0aCwiaW5wdXRzL0xvRGVuc19nZW9kaXN0X3BhLnJkcyIpKSkNCmBgYA0KDQojIyMjIE9yZGluYXRpb24NCg0KIyMjIyBtb25vTURTIC0gZG9uJ3QgcnVuIGlmIGxvYWRpbmcgc2F2ZWQgb2JqZWN0DQpgYGB7cn0NCg0KU3lzLnRpbWUoKQ0KZCA8LSAyDQptZHNfcGEgPC0gbGlzdCgpDQoNCmZvciAoaSBpbiAxOjEwMCkgew0KICBtZHNfcGFbW2ldXTwtbW9ub01EUyhnZW9kaXN0X3BhLA0KICAgICAgICAgICAgICAgICAgICBtYXRyaXgoYyhydW5pZihkaW0ob3R1X3BhKVsxXSpkKSksDQogICAgICAgICAgICAgICAgICAgICAgICAgICBucm93ID0gZGltKG90dV9wYSlbMV0pLA0KICAgICAgICAgICAgICAgICAgICBrID0gZCwNCiAgICAgICAgICAgICAgICAgICAgbW9kZWwgPSAiZ2xvYmFsIiwNCiAgICAgICAgICAgICAgICAgICAgbWF4aXQgPSAyMDAwLA0KICAgICAgICAgICAgICAgICAgICBzbWluID0gMWUtNywNCiAgICAgICAgICAgICAgICAgICAgc2Zncm1pbiA9IDFlLTcpDQp9DQoNCmBgYA0KDQojIyMjIyBTYXZlIHRoZSByZXN1bHQgLSBkb24ndCBydW4gaWYgbG9hZGluZyBzYXZlZCBvYmplY3QNCmBgYHtyfQ0Kc2F2ZVJEUyhtZHNfcGEsDQogICAgICAgIGZpbGUgPSBmaWxlLnBhdGgoZGF0YVBhdGgsImlucHV0cy9Mb0RlbnNfbWRzX3BhLnJkcyIpKQ0KYGBgDQoNCg0KIyMjIyMgQmVzdCBubWRzIHNvbHV0aW9uIC0gUEENCk1ha2UgZnVuY3Rpb24/DQpgYGB7ciwgZXZhbD1UUlVFfQ0KIyBMb2FkaW5nIGdlb2Rpc3Qgb2JqZWN0DQojIGdlb2Rpc3RfbmZpIDwtIHJlYWRSRFMoZmlsZSA9ICIuLi9TcGxpdFJldjJfZ2VvZGlzdF9wYV9mdWxsLnJkcyIpDQpgYGANCg0KDQpgYGB7cn0NCiMgTG9hZGluZyBtZHMgcmVzdWx0cw0KIyBtZHMgPC0gcmVhZFJEUyhmaWxlID0gIi4uL1NwbGl0UmV2Ml9tZHNfcGFfZnVsbC5yZHMiKQ0KDQojIyBFeHRyYWN0aW5nIHRoZSBzdHJlc3Mgb2YgZWFjaCBubWRzIGl0ZXJhdGlvbg0KbWRzX3N0cmVzc19wYTwtdW5saXN0KGxhcHBseShtZHNfcGEsIGZ1bmN0aW9uKHYpe3ZbWzIyXV19KSkgDQoNCm9yZGVyZWRfcGEgPC1vcmRlcihtZHNfc3RyZXNzX3BhKQ0KDQojIyBCZXN0LCBzZWNvbmQgYmVzdCwgYW5kIHdvcnN0IHNvbHV0aW9uDQptZHNfc3RyZXNzX3BhW29yZGVyZWRfcGFbMV1dDQptZHNfc3RyZXNzX3BhW29yZGVyZWRfcGFbMl1dDQptZHNfc3RyZXNzX3BhW29yZGVyZWRfcGFbMTBdXQ0KDQojIyBTY2FsaW5nIG9mIGF4ZXMgdG8gaGFsZiBjaGFuZ2UgdW5pdHMgYW5kIHZhcmltYXggcm90YXRpb24gYnkgcG9zdE1EUw0KbWRzX2Jlc3RfcGE8LXBvc3RNRFMobWRzX3BhW1tvcmRlcmVkX3BhWzFdXV0sDQogICAgICAgICAgICAgICAgICBnZW9kaXN0X3BhLCANCiAgICAgICAgICAgICAgICAgIHBjID0gVFJVRSwgDQogICAgICAgICAgICAgICAgICBoYWxmY2hhbmdlID0gVFJVRSwgDQogICAgICAgICAgICAgICAgICB0aHJlc2hvbGQgPSBlcCkgICAgICMgSXMgdGhpcyB0aHJlc2hvbGQgcmVsYXRlZCB0byB0aGUgZXBzaWxvbiBhYm92ZT8NCm1kc19iZXN0X3BhDQoNCm1kc19zZWNiZXN0X3BhIDwtIHBvc3RNRFMobWRzX3BhW1tvcmRlcmVkX3BhWzJdXV0sDQogICAgICAgICAgICAgICAgICAgICAgICAgIGdlb2Rpc3RfcGEsIA0KICAgICAgICAgICAgICAgICAgICAgICAgICBwYyA9IFRSVUUsIA0KICAgICAgICAgICAgICAgICAgICAgICAgICBoYWxmY2hhbmdlID0gVFJVRSwgDQogICAgICAgICAgICAgICAgICAgICAgICAgIHRocmVzaG9sZCA9IGVwKQ0KbWRzX3NlY2Jlc3RfcGENCg0KIyMgUHJvY3J1c3RlcyBjb21wYXJpc29ucw0KcHJvY3JfcGEgPC0gcHJvY3J1c3RlcyhtZHNfYmVzdF9wYSwNCiAgICAgICAgICAgICAgICAgICAgbWRzX3NlY2Jlc3RfcGEsDQogICAgICAgICAgICAgICAgICAgIHBlcm11dGF0aW9ucz05OTkpDQpwcm90ZXN0KG1kc19iZXN0X3BhLA0KICAgICAgICBtZHNfc2VjYmVzdF9wYSwNCiAgICAgICAgcGVybXV0YXRpb25zPTk5OSkNCg0KcGxvdChwcm9jcl9wYSkNCg0KcG5nKGZpbGU9ZmlsZS5wYXRoKGRhdGFQYXRoLCJvdXRwdXRzL0xvRGVuc19wcm9jcnVzdGVzX3BhLnBuZyIpLCB3aWR0aD0xMDAwLCBoZWlnaHQ9NzAwKQ0KcGxvdChwcm9jcl9wYSkNCmRldi5vZmYoKQ0KYGBgDQoNCiMjIyMjIENvcnJlbGF0aW9uIG9mIGF4aXM6IERDQSB2cyBOTURTIC0gUEENCmBgYHtyfQ0KIyBFeHRyYWN0aW5nIG9yZGluYXRpb24gYXhpcw0KYXggPC0gMg0KYXhpc19wYSA8LSBjYmluZChtZHNfYmVzdF9wYSRwb2ludHMsDQogICAgICAgICAgICAgICAgIHNjb3JlcyhkY2FfcGEsDQogICAgICAgICAgICAgICAgICAgICAgICBkaXNwbGF5ID0gInNpdGVzIiwNCiAgICAgICAgICAgICAgICAgICAgICAgIG9yaWdpbiA9IFRSVUUpWywgMTpheF0pDQoNCmdnY29ycihheGlzX3BhLCANCiAgICAgICBtZXRob2Q9YygiZXZlcnl0aGluZyIsImtlbmRhbGwiKSwgDQogICAgICAgbGFiZWwgPSBUUlVFLA0KICAgICAgIGxhYmVsX3NpemUgPSAzLCANCiAgICAgICBsYWJlbF9jb2xvciA9ICJibGFjayIsICANCiAgICAgICBuYnJlYWtzID0gOCwNCiAgICAgICBsYWJlbF9yb3VuZCA9IDMsDQogICAgICAgbG93ID0gInJlZCIsDQogICAgICAgbWlkID0gIndoaXRlIiwNCiAgICAgICBoaWdoID0gImdyZWVuIikNCg0KDQpgYGANCg0KIyMjIyMgU2F2ZSB0aGUgZmlndXJlDQpgYGB7cn0NCmdnc2F2ZShmaWxlbmFtZSA9IGZpbGUucGF0aChkYXRhUGF0aCwib3V0cHV0cy9Mb0RlbnNfY29ycmVsYXRpb25QQ0F2c05NRFNfUEEucG5nIiksDQogICAgICAgZGV2aWNlID0gInBuZyIsDQogICAgICAgZHBpPTMwMCApDQoNCiMgU3dpdGNoaW5nIGRpcmVjdGlvbiBvZiBOTURTMQ0KIyBtZHNfYmVzdCRwb2ludHNbLCAxXSA8LSAtbWRzX2Jlc3QkcG9pbnRzWywgMV0NCmBgYA0KDQojIyMgRENBIFI2DQpgYGB7cn0NCmRjYV9yNiA8LSBkZWNvcmFuYSh2ZWcgPSBvdHVfNikNCg0KcHJpbnQoZGNhX3I2LCBoZWFkPVQpDQpgYGANCg0KDQojIyMgR05NRFMgUjYNCiMjIyMgRGlzdGFuY2VzIC0gZG9uJ3QgcnVuIGlmIGxvYWRpbmcgc2F2ZWQgb2JqZWN0DQpgYGB7cn0NCiMjIEJyYXktQ3VydGlzDQpkaXN0X3I2IDwtIHZlZ2Rpc3QoeCA9IG90dV82LCBtZXRob2QgPSAiYnJheSIpDQoNCiMjIEdlb2Rpc3QNCmVwIDwtIDAuODAgICAgICMgZXBzaWxvbg0KZ2VvZGlzdF9yNiA8LSBpc29tYXBkaXN0KGRpc3QgPSBkaXN0X3I2LCBlcHNpbG9uID0gZXApDQpgYGANCg0KIyMjIyMgU2F2ZSB0aGUgcmVzdWx0IC0gZG9uJ3QgcnVuIGlmIGxvYWRpbmcgc2F2ZWQgb2JqZWN0DQpgYGB7cn0NCnNhdmVSRFMoZ2VvZGlzdF9yNiwNCiAgICAgICAgZmlsZSA9IChmaWxlLnBhdGgoZGF0YVBhdGgsImlucHV0cy9Mb0RlbnNfZ2VvZGlzdF9yNi5yZHMiKSkpDQpgYGANCg0KDQojIyMjIE9yZGluYXRpb24gDQoNCiMjIyMjIDIwMCByZXBzIC0gZG9uJ3QgcnVuIGlmIGxvYWRpbmcgc2F2ZWQgb2JqZWN0DQp0b29rIDIwbWlucw0KYGBge3J9DQojIG1vbm9NRFMNCmQgPC0gMg0KbWRzX3I2IDwtIGxpc3QoKQ0KDQpTeXMudGltZSgpDQpmb3IgKGkgaW4gMToyMDApIHsNCiAgbWRzX3I2W1tpXV08LW1vbm9NRFMoZ2VvZGlzdF9yNiwNCiAgICAgICAgICAgICAgICAgICAgbWF0cml4KGMocnVuaWYoZGltKG90dV82KVsxXSpkKSksDQogICAgICAgICAgICAgICAgICAgICAgICAgICBucm93ID0gZGltKG90dV82KVsxXSksDQogICAgICAgICAgICAgICAgICAgIGsgPSBkLA0KICAgICAgICAgICAgICAgICAgICBtb2RlbCA9ICJnbG9iYWwiLA0KICAgICAgICAgICAgICAgICAgICBtYXhpdCA9IDIwMDAsDQogICAgICAgICAgICAgICAgICAgIHNtaW4gPSAxZS03LA0KICAgICAgICAgICAgICAgICAgICBzZmdybWluID0gMWUtNykNCn0NClN5cy50aW1lKCkNCmBgYA0KDQojIyMjIyBTYXZlIHRoZSByZXN1bHQgLSBkb24ndCBydW4gaWYgbG9hZGluZyBzYXZlZCBvYmplY3QNCmBgYHtyfQ0Kc2F2ZVJEUyhtZHNfcjYsDQogICAgICAgIGZpbGUgPSBmaWxlLnBhdGgoZGF0YVBhdGgsImlucHV0cy9Mb0RlbnNfbWRzX3I2LnJkcyIpKSANCmBgYA0KDQojIyMjIyBCZXN0IG5tZHMgc29sdXRpb24gcjYgMjAwIHJlcCANCmBgYHtyLCBldmFsPVRSVUV9DQojIExvYWRpbmcgZ2VvZGlzdCBvYmplY3QNCiMgZ2VvZGlzdF9uZmkgPC0gcmVhZFJEUyhmaWxlID0gIi4uL1NwbGl0UmV2Ml9nZW9kaXN0X25maS5yZHMiKQ0KDQojIExvYWRpbmcgbWRzIHJlc3VsdHMNCiMgbWRzIDwtIHJlYWRSRFMoZmlsZSA9ICIuLi9TcGxpdFJldjJfbWRzLnJkcyIpDQoNCiMjIEV4dHJhY3RpbmcgdGhlIHN0cmVzcyBvZiBlYWNoIG5tZHMgaXRlcmF0aW9uDQptZHNfc3RyZXNzX3I2PC11bmxpc3QobGFwcGx5KG1kc19yNiwgZnVuY3Rpb24odil7dltbMjJdXX0pKSANCg0Kb3JkZXJlZF9yNiA8LW9yZGVyKG1kc19zdHJlc3NfcjYpDQoNCiMjIEJlc3QsIHNlY29uZCBiZXN0LCBhbmQgd29yc3Qgc29sdXRpb24NCm1kc19zdHJlc3NfcjZbb3JkZXJlZF9yNlsxXV0NCm1kc19zdHJlc3NfcjZbb3JkZXJlZF9yNlsyXV0NCm1kc19zdHJlc3NfcjZbb3JkZXJlZF9yNlsxMDBdXQ0KDQojIyBTY2FsaW5nIG9mIGF4ZXMgdG8gaGFsZiBjaGFuZ2UgdW5pdHMgYW5kIHZhcmltYXggcm90YXRpb24gYnkgcG9zdE1EUw0KbWRzX2Jlc3RfcjY8LXBvc3RNRFMobWRzX3I2W1tvcmRlcmVkX3I2WzFdXV0sDQogICAgICAgICAgICAgICAgICAgICBnZW9kaXN0X3I2LCANCiAgICAgICAgICAgICAgICAgICAgIHBjID0gVFJVRSwgDQogICAgICAgICAgICAgICAgICAgICBoYWxmY2hhbmdlID0gVFJVRSwgDQogICAgICAgICAgICAgICAgICAgICB0aHJlc2hvbGQgPSBlcCkgICAgICMgSXMgdGhpcyB0aHJlc2hvbGQgcmVsYXRlZCB0byB0aGUgZXBzaWxvbiBhYm92ZT8NCm1kc19iZXN0X3I2DQoNCm1kc19zZWNiZXN0X3I2PC1wb3N0TURTKG1kc19yNltbb3JkZXJlZF9yNlsyXV1dLA0KICAgICAgICAgICAgICAgICAgICAgICAgZ2VvZGlzdF9yNiwgDQogICAgICAgICAgICAgICAgICAgICAgICBwYyA9IFRSVUUsIA0KICAgICAgICAgICAgICAgICAgICAgICAgaGFsZmNoYW5nZSA9IFRSVUUsIA0KICAgICAgICAgICAgICAgICAgICAgICAgdGhyZXNob2xkID0gZXApDQptZHNfc2VjYmVzdF9yNg0KDQojIyBQcm9jcnVzdGVzIGNvbXBhcmlzb25zDQpwcm9jcl9yNiA8LSBwcm9jcnVzdGVzKG1kc19iZXN0X3I2LA0KICAgICAgICAgICAgICAgICAgICAgICBtZHNfc2VjYmVzdF9yNiwNCiAgICAgICAgICAgICAgICAgICAgICAgcGVybXV0YXRpb25zPTk5OSkNCnByb3Rlc3QobWRzX2Jlc3RfcjYsDQogICAgICAgIG1kc19zZWNiZXN0X3I2LA0KICAgICAgICBwZXJtdXRhdGlvbnM9OTk5KQ0KDQpwbG90KHByb2NyX3I2KQ0KDQpwbmcoZmlsZS5wYXRoKGRhdGFQYXRoLCJvdXRwdXRzL0xvRGVuc19wcm9jcnVzdGVzX3I2LnBuZyIpLCB3aWR0aD0xMDAwLCBoZWlnaHQ9NzAwLCkgI2FkZGVkIDEwMDANCnBsb3QocHJvY3JfcjYpDQpkZXYub2ZmKCkNCmBgYA0KDQogIyMjIyAxMDAwIHJlcHMgLSBkb24ndCBydW4gaWYgbG9hZGluZyBzYXZlZCBvYmplY3QgDQogQ3VycmVudGx5IGNvbW1lbnRlZCBvdXQgYXMgaXQgZ2FpbmVkIG5vdGhpbmcgYnV0IHRvb2sgZXh0cmEgdGltZS4gQ2FuIGJlIHJlbW92ZWQgaW4gZHVlIGNvdXJzZS4NCjwhLS0gdG9vayAxaHIgLS0+DQo8IS0tIGBgYHtyfSAtLT4NCjwhLS0gIyBtb25vTURTIC0tPg0KPCEtLSBkIDwtIDIgLS0+DQo8IS0tIG1kc19yNl8xMDAwIDwtIGxpc3QoKSNvYmogYWRkZWRfMTAwMCAtLT4NCg0KPCEtLSBTeXMudGltZSgpIC0tPg0KPCEtLSBmb3IgKGkgaW4gMToxMDAwKSB7I2NoYW5nZSBmcm9tIDIwMCAoMjAwIHJhIGluIDIwbWluLCAxMDAwIHJhbiBpbiAxaHIpIC0tPg0KPCEtLSAgIG1kc19yNl8xMDAwW1tpXV08LW1vbm9NRFMoZ2VvZGlzdF9yNiwjb2JqIGFkZGVkIF8xMDAwIC0tPg0KPCEtLSAgICAgICAgICAgICAgICAgICAgIG1hdHJpeChjKHJ1bmlmKGRpbShvdHVfNilbMV0qZCkpLCAtLT4NCjwhLS0gICAgICAgICAgICAgICAgICAgICAgICAgICAgbnJvdyA9IGRpbShvdHVfNilbMV0pLCAtLT4NCjwhLS0gICAgICAgICAgICAgICAgICAgICBrID0gZCwgLS0+DQo8IS0tICAgICAgICAgICAgICAgICAgICAgbW9kZWwgPSAiZ2xvYmFsIiwgLS0+DQo8IS0tICAgICAgICAgICAgICAgICAgICAgbWF4aXQgPSAyMDAwLCAtLT4NCjwhLS0gICAgICAgICAgICAgICAgICAgICBzbWluID0gMWUtNywgLS0+DQo8IS0tICAgICAgICAgICAgICAgICAgICAgc2Zncm1pbiA9IDFlLTcpIC0tPg0KPCEtLSB9IC0tPg0KPCEtLSBTeXMudGltZSgpIC0tPg0KPCEtLSBgYGAgLS0+DQoNCjwhLS0gIyMjIyMgU2F2ZSB0aGUgcmVzdWx0IC0gZG9uJ3QgcnVuIGlmIGxvYWRpbmcgc2F2ZWQgb2JqZWN0IC0tPg0KPCEtLSBgYGB7cn0gLS0+DQo8IS0tIHNhdmVSRFMobWRzX3I2XzEwMDAsICNvYmogYWRkZWRfMTAwMCAtLT4NCjwhLS0gICAgICAgICBmaWxlID0gZmlsZS5wYXRoKGRhdGFQYXRoLCJpbnB1dHMvbWRzX3I2XzEwMDByZXAucmRzIikpICNjaGFuZ2UgZnJvbSAyMDAgLS0+DQo8IS0tIGBgYCAtLT4NCg0KPCEtLSAjIyMjIyBCZXN0IG5tZHMgc29sdXRpb24gcjYgMTAwMCByZXAgLS0+DQo8IS0tIGBgYHtyLCBldmFsPVRSVUV9IC0tPg0KPCEtLSAjIExvYWRpbmcgZ2VvZGlzdCBvYmplY3QgLS0+DQo8IS0tICMgZ2VvZGlzdF9uZmkgPC0gcmVhZFJEUyhmaWxlID0gIi4uL1NwbGl0UmV2Ml9nZW9kaXN0X25maS5yZHMiKSAtLT4NCg0KPCEtLSAjIExvYWRpbmcgbWRzIHJlc3VsdHMgLS0+DQo8IS0tICMgbWRzIDwtIHJlYWRSRFMoZmlsZSA9ICIuLi9TcGxpdFJldjJfbWRzLnJkcyIpIC0tPg0KDQo8IS0tICMjIEV4dHJhY3RpbmcgdGhlIHN0cmVzcyBvZiBlYWNoIG5tZHMgaXRlcmF0aW9uIC0tPg0KPCEtLSBtZHNfc3RyZXNzX3I2XzEwMDA8LXVubGlzdChsYXBwbHkobWRzX3I2XzEwMDAsIGZ1bmN0aW9uKHYpe3ZbWzIyXV19KSkgI2FkZWVkIDEwMDAgLS0+DQoNCjwhLS0gb3JkZXJlZF9yNl8xMDAwIDwtb3JkZXIobWRzX3N0cmVzc19yNl8xMDAwKSAtLT4NCg0KPCEtLSAjIyBCZXN0LCBzZWNvbmQgYmVzdCwgYW5kIHdvcnN0IHNvbHV0aW9uIC0tPg0KPCEtLSBtZHNfc3RyZXNzX3I2XzEwMDBbb3JkZXJlZF9yNl8xMDAwWzFdXSAtLT4NCjwhLS0gbWRzX3N0cmVzc19yNl8xMDAwW29yZGVyZWRfcjZfMTAwMFsyXV0gLS0+DQo8IS0tIG1kc19zdHJlc3NfcjZfMTAwMFtvcmRlcmVkX3I2XzEwMDBbMTAwXV0gLS0+DQoNCjwhLS0gIyMgU2NhbGluZyBvZiBheGVzIHRvIGhhbGYgY2hhbmdlIHVuaXRzIGFuZCB2YXJpbWF4IHJvdGF0aW9uIGJ5IHBvc3RNRFMgLS0+DQo8IS0tIG1kc19iZXN0X3I2XzEwMDA8LXBvc3RNRFMobWRzX3I2XzEwMDBbW29yZGVyZWRfcjZfMTAwMFsxXV1dLCAtLT4NCjwhLS0gICAgICAgICAgICAgICAgICAgICAgZ2VvZGlzdF9yNiwgIC0tPg0KPCEtLSAgICAgICAgICAgICAgICAgICAgICBwYyA9IFRSVUUsICAtLT4NCjwhLS0gICAgICAgICAgICAgICAgICAgICAgaGFsZmNoYW5nZSA9IFRSVUUsICAtLT4NCjwhLS0gICAgICAgICAgICAgICAgICAgICAgdGhyZXNob2xkID0gZXApICAgICAjIElzIHRoaXMgdGhyZXNob2xkIHJlbGF0ZWQgdG8gdGhlIGVwc2lsb24gYWJvdmU/IC0tPg0KPCEtLSBtZHNfYmVzdF9yNl8xMDAwIC0tPg0KDQo8IS0tIG1kc19zZWNiZXN0X3I2XzEwMDA8LXBvc3RNRFMobWRzX3I2XzEwMDBbW29yZGVyZWRfcjZfMTAwMFsyXV1dLCAtLT4NCjwhLS0gICAgICAgICAgICAgICAgICAgICAgICAgZ2VvZGlzdF9yNiwgIC0tPg0KPCEtLSAgICAgICAgICAgICAgICAgICAgICAgICBwYyA9IFRSVUUsICAtLT4NCjwhLS0gICAgICAgICAgICAgICAgICAgICAgICAgaGFsZmNoYW5nZSA9IFRSVUUsICAtLT4NCjwhLS0gICAgICAgICAgICAgICAgICAgICAgICAgdGhyZXNob2xkID0gZXApIC0tPg0KPCEtLSBtZHNfc2VjYmVzdF9yNl8xMDAwIC0tPg0KDQo8IS0tICMjIFByb2NydXN0ZXMgY29tcGFyaXNvbnMgLS0+DQo8IS0tIHByb2NyX3I2XzEwMDAgPC0gcHJvY3J1c3RlcyhtZHNfYmVzdF9yNl8xMDAwLCAtLT4NCjwhLS0gICAgICAgICAgICAgICAgICAgICAgICBtZHNfc2VjYmVzdF9yNl8xMDAwLCAtLT4NCjwhLS0gICAgICAgICAgICAgICAgICAgICAgICBwZXJtdXRhdGlvbnM9OTk5KSAtLT4NCjwhLS0gcHJvdGVzdChtZHNfYmVzdF9yNl8xMDAwLCAtLT4NCjwhLS0gICAgICAgICBtZHNfc2VjYmVzdF9yNl8xMDAwLCAtLT4NCjwhLS0gICAgICAgICBwZXJtdXRhdGlvbnM9OTk5KSAtLT4NCg0KPCEtLSBwbG90KHByb2NyX3I2XzEwMDApIC0tPg0KDQo8IS0tIHBuZyhmaWxlLnBhdGgoZGF0YVBhdGgsIm91dHB1dHMvcHJvY3J1c3Rlc19yNl8xMDAwLnBuZyIpLCB3aWR0aD0xMDAwLCBoZWlnaHQ9NzAwLCkgI2FkZGVkIDEwMDAgLS0+DQo8IS0tIHBsb3QocHJvY3JfcjZfMTAwMCkgLS0+DQo8IS0tIGRldi5vZmYoKSAtLT4NCjwhLS0gYGBgIC0tPg0KDQoNCg0KIyMjIyMgQ29ycmVsYXRpb24gb2YgYXhpczogRENBIHZzIE5NRFMgDQpyZXRhaW4gdGhlIDIwMCByZXAgdmVyc2lvbg0KYGBge3J9DQojIEV4dHJhY3Rpbmcgb3JkaW5hdGlvbiBheGlzDQpheCA8LSAyDQpheGlzX3I2IDwtIGNiaW5kKG1kc19iZXN0X3I2JHBvaW50cywNCiAgICAgICAgICAgICAgICAgc2NvcmVzKGRjYV9yNiwNCiAgICAgICAgICAgICAgICAgICAgICAgIGRpc3BsYXkgPSAic2l0ZXMiLA0KICAgICAgICAgICAgICAgICAgICAgICAgb3JpZ2luID0gVFJVRSlbLCAxOmF4XSkNCg0KZ2djb3JyKGF4aXNfcjYsIA0KICAgICAgIG1ldGhvZD1jKCJldmVyeXRoaW5nIiwia2VuZGFsbCIpLCANCiAgICAgICBsYWJlbCA9IFRSVUUsDQogICAgICAgbGFiZWxfc2l6ZSA9IDMsIA0KICAgICAgIGxhYmVsX2NvbG9yID0gImJsYWNrIiwgIA0KICAgICAgIG5icmVha3MgPSA4LA0KICAgICAgIGxhYmVsX3JvdW5kID0gMywNCiAgICAgICBsb3cgPSAicmVkIiwNCiAgICAgICBtaWQgPSAid2hpdGUiLA0KICAgICAgIGhpZ2ggPSAiZ3JlZW4iKQ0KDQpgYGANCg0KIyMjIyMgU2F2ZSB0aGUgZmlndXJlDQpgYGB7cn0NCmdnc2F2ZShmaWxlbmFtZSA9IGZpbGUucGF0aChkYXRhUGF0aCwib3V0cHV0cy9Mb0RlbnNfY29ycmVsYXRpb25QQ0F2c05NRFNfcjYucG5nIiksDQogICAgICAgZGV2aWNlID0gInBuZyIsDQogICAgICAgZHBpPTMwMCApDQoNCiMgU3dpdGNoaW5nIGRpcmVjdGlvbiBvZiBOTURTMQ0KIyBtZHNfYmVzdCRwb2ludHNbLCAxXSA8LSAtbWRzX2Jlc3QkcG9pbnRzWywgMV0NCmBgYA0KDQojIyMgUGxvdHRpbmcgRENBICYgR05NRFMNCmBgYHtyfQ0KIyMgQWRkaW5nIHNjb3JlcyB0byBkYXRhIGZyYW1lDQpvdHVfNiRnbm1kczEgPC0gbWRzX2Jlc3RfcjYkcG9pbnRzWywgMV0NCm90dV82JGdubWRzMiA8LSBtZHNfYmVzdF9yNiRwb2ludHNbLCAyXQ0Kb3R1XzYkZGNhMSA8LSBzY29yZXMoZGNhX3I2LCBkaXNwbGF5ID0gInNpdGVzIiwgb3JpZ2luID0gVFJVRSlbLCAxXQ0Kb3R1XzYkZGNhMiA8LSBzY29yZXMoZGNhX3I2LCBkaXNwbGF5ID0gInNpdGVzIiwgb3JpZ2luID0gVFJVRSlbLCAyXQ0KDQpwX2dubWRzX3I2IDwtIGdncGxvdChkYXRhID0gb3R1XzYsDQogICAgICAgICAgICAgICAgICAgICBhZXMoeCA9IGdubWRzMSwNCiAgICAgICAgICAgICAgICAgICAgICAgICB5ID0gZ25tZHMyKSkgKw0KICB0aGVtZV9jbGFzc2ljKCkgKw0KICBjb29yZF9maXhlZCgpICsNCiAgZ2d0aXRsZSgiR05NRFMiLA0KICAgICAgICAgIHN1YnRpdGxlID0gIkZpcnN0IHJ1biIpICsNCiAgZ2VvbV9wb2ludChjb2xvdXIgPSAicmVkIikgKw0KICBnZW9tX3ZsaW5lKHhpbnRlcmNlcHQgPSAwLA0KICAgICAgICAgICAgIGxpbmV0eXBlID0gMiwNCiAgICAgICAgICAgICBjb2xvdXIgPSAibGlnaHRncmF5IikgKw0KICBnZW9tX2hsaW5lKHlpbnRlcmNlcHQgPSAwLA0KICAgICAgICAgICAgIGxpbmV0eXBlID0gMiwNCiAgICAgICAgICAgICBjb2xvdXIgPSAibGlnaHRncmF5IikNCg0KcF9kY2FfcjYgPC0gZ2dwbG90KGRhdGEgPSBvdHVfNiwNCiAgICAgICAgICAgICAgICAgICBhZXMoeCA9IGRjYTEsDQogICAgICAgICAgICAgICAgICAgICAgIHkgPSBkY2EyKSkgKw0KICB0aGVtZV9jbGFzc2ljKCkgKw0KICBjb29yZF9maXhlZCgpICsNCiAgZ2d0aXRsZSgiRENBIiwNCiAgICAgICAgICBzdWJ0aXRsZSA9ICJGaXJzdCBydW4iKSArDQogIGdlb21fcG9pbnQoY29sb3VyID0gInJlZCIpICsNCiAgZ2VvbV92bGluZSh4aW50ZXJjZXB0ID0gMCwNCiAgICAgICAgICAgICBsaW5ldHlwZSA9IDIsDQogICAgICAgICAgICAgY29sb3VyID0gImxpZ2h0Z3JheSIpICsNCiAgZ2VvbV9obGluZSh5aW50ZXJjZXB0ID0gMCwNCiAgICAgICAgICAgICBsaW5ldHlwZSA9IDIsDQogICAgICAgICAgICAgY29sb3VyID0gImxpZ2h0Z3JheSIpDQoNCnBfZ25tZHNfcjYgKyBwX2RjYV9yNg0KDQoNCmBgYA0KDQojIyMjIyBTYXZlIHRoZSBmaWd1cmUNCmBgYHtyfQ0KZ2dzYXZlKGZpbGUucGF0aChkYXRhUGF0aCwib3V0cHV0cy9Mb0RlbnNfZ25tZHNfZGNhX3I2LnBuZyIpLA0KICAgICAgIGRldmljZSA9ICJwbmciLCANCiAgICAgICBkcGk9MzAwKQ0KYGBgDQoNCiMjIFNwZWNpZXMtZW52aXJvbm1lbnQgcmVsYXRpb25zaGlwcw0KIyMjIFNlbGVjdGluZyBvcmRpbmF0aW9uDQpgYGB7cn0NCm9yZCA8LSBtZHNfYmVzdF9yNg0KDQojIyBBeGlzIHNjb3JlcyBpZiBzZWxlY3RlZCBvcmQgaXMgR05NRFMNCmF4aXMgPC0gb3JkJHBvaW50cyAlPiUgYXMuZGF0YS5mcmFtZQ0KDQojIyBBeGlzIHNjb3JlcyBpZiBzZWxlY3RlZCBvcmQgaXMgRENBDQojIGF4aXMgPC0gc2NvcmVzKG9yZCwNCiMgICAgICAgICAgICAgICAgZGlzcGxheSA9ICJzaXRlcyIsDQojICAgICAgICAgICAgICAgIG9yaWdpbiA9IFRSVUUpWywgMTpheF0pDQoNCmBgYA0KDQojIyMgQ3JlYXRlIGFkZGl0aW9uYWwgdmFyaWFibGVzDQpkZWNpZGVkIHRvIG1ha2UgTUxELWJhdGh5IHZhcnMNCmBgYHtyfQ0KZW52PC1lbnYgJT4lIA0KICBtdXRhdGUgKCJNTERtZWFuX2JhdGh5Ij1NTERtZWFuX1JvYmluc29uLShiYXRoeSotMSksDQogICAgICAgICAgIk1MRG1pbl9iYXRoeSI9TUxEbWluX1JvYmluc29uLShiYXRoeSotMSksDQogICAgICAgICAgIk1MRG1heF9iYXRoeSI9TUxEbWF4X1JvYmluc29uLShiYXRoeSotMSkpDQoNCmVudiRNTERtZWFuX2JhdGh5PC1jdXQoZW52JE1MRG1lYW5fYmF0aHksIA0KICAgICAgYnJlYWtzPWMoLTI1NjAsIC0yMCwyMCwxMzApLCNjaGVja2VkIHJhbmdlIG9mIHZhbHVlcyBmaXJzdCAobWluIC0yNTU0LCBtYXggMTIzKQ0KICAgICAgbGFiZWxzPWMoJ2JlbG93TUxEJywnb25QeWNubycsJ2luTWl4TGF5ZXInKSkNCmVudiRNTERtaW5fYmF0aHk8LWN1dChlbnYkTUxEbWluX2JhdGh5LCANCiAgICAgIGJyZWFrcz1jKC0yNTYwLCAtMjAsMjAsMTMwKSwjY2hlY2tlZCByYW5nZSBvZiB2YWx1ZXMgZmlyc3QgKG1pbiAtMjU1NCwgbWF4IDEyMykNCiAgICAgIGxhYmVscz1jKCdiZWxvd01MRCcsJ29uUHljbm8nLCdpbk1peExheWVyJykpDQplbnYkTUxEbWF4X2JhdGh5PC1jdXQoZW52JE1MRG1heF9iYXRoeSwgDQogICAgICBicmVha3M9YygtMjU2MCwgLTIwLDIwLDEzMCksI2NoZWNrZWQgcmFuZ2Ugb2YgdmFsdWVzIGZpcnN0IChtaW4gLTI1NTQsIG1heCAxMjMpDQogICAgICBsYWJlbHM9YygnYmVsb3dNTEQnLCdvblB5Y25vJywnaW5NaXhMYXllcicpKQ0KDQplbnYkc3dEZW5zUm9iX2F2czwtc3dSaG8oc2FsaW5pdHk9ZW52JFNtZWFuX1JvYmluc29uLA0KICAgICAgICAgICAgICAgICAgICAgICAgIHRlbXBlcmF0dXJlPWVudiRUbWVhbl9Sb2JpbnNvbiwNCiAgICAgICAgICAgICAgICAgICAgICAgICBwcmVzc3VyZT0oZW52JGJhdGh5Ki0xKSwNCiAgICAgICAgICAgICAgICAgICAgICAgICBlb3M9InVuZXNjbyIpDQoNCmBgYA0KDQoNCg0KIyMjIENvcnJlbGF0aW9uIG9yZGluYXRpb24gYXhlcyBhbmQgZW52aXJvbm1lbnRhbCB2YXJpYWJsZXMNCiMjIyMgUmVtb3Zpbmcgbm9uLWVudiB2YXJzDQpgYGB7cn0NCmVudl9jb250PC1lbnYlPiUgc2VsZWN0KC1jKGxhbmRzY2FwZSxzZWRjbGFzcyxnbW9ycGgsIE1MRG1lYW5fYmF0aHksIE1MRG1heF9iYXRoeSwgTUxEbWluX2JhdGh5LCAjY2F0ZWdvcmljYWwNCiAgICAgICAgICAgICAgICAgICAgICAgICAgIG9wdGlvbmFsLCAjbm90IGEgdmFyDQogICAgICAgICAgICAgICAgICAgICAgICAgICBNTERtYXhfUm9iaW5zb24sIE1MRG1lYW5fUm9iaW5zb24sIE1MRG1pbl9Sb2JpbnNvbiwgICNyZXBsYWNlZCBieSBuZXcgdmFycw0KICAgICAgICAgICAgICAgICAgICAgICAgICAgTUxEc2RfUm9iaW5zb24gI25vdCBtZWFuaW5nZnVsDQogICAgICAgICAgICAgICAgICAgICAgICAgICApKQ0KZW52X2NvbnQ8LWVudl9jb250JT4lIG11dGF0ZV9pZihpcy5pbnRlZ2VyLGFzLm51bWVyaWMpDQplbnZfY29yciA8LSBlbnZfY29udCAlPiUgc2VsZWN0KC1jKFNhbXBJRCkpDQoNCiMgZW52X2NvcnIkY29vcmRzLngxPC1hcy5udW1lcmljKGVudl9jb3JyJGNvb3Jkcy54MSkNCiMgZW52X2NvcnIkY29vcmRzLngyPC1hcy5udW1lcmljKGVudl9jb3JyJGNvb3Jkcy54MikNCg0KZW52X2NvcnJbKCFpcy5udW1lcmljKGVudl9jb3JyKSksXQ0KYGBgDQoNCiMjIyMgQ29ycmVsYXRpb25zDQpgYGB7cn0NCiMgVmVjdG9yIHRvIGhvbGQgY29ycmVsYXRpb25zDQpjb3JfYXgxIDwtIE5VTEwNCmNvcl9heDIgPC0gTlVMTA0KcHZfYXgxIDwtIE5VTEwNCnB2X2F4MiA8LSBOVUxMDQoNCiMgTk1EUzENCmZvciggaSBpbiBzZXEobGVuZ3RoKGVudl9jb3JyKSkpIHsNCiAgY3QuaSA8LSBjb3IudGVzdChheGlzJE1EUzEsDQogICAgICAgICAgICAgICAgICAgZW52X2NvcnJbLCBpXSwNCiAgICAgICAgICAgICAgICAgICBtZXRob2QgPSAia2VuZGFsbCIpDQogIGNvcl9heDFbaV0gPC0gY3QuaSRlc3RpbWF0ZQ0KICBwdl9heDFbaV0gPC0gY3QuaSRwLnZhbHVlDQp9DQoNCiMgTk1EUzINCmZvciggaSBpbiBzZXEobGVuZ3RoKGVudl9jb3JyKSkpIHsNCiAgY3QuaSA8LSBjb3IudGVzdChheGlzJE1EUzIsDQogICAgICAgICAgICAgICAgICAgZW52X2NvcnJbLCBpXSwNCiAgICAgICAgICAgICAgICAgICBtZXRob2QgPSAia2VuZGFsbCIpDQogIGNvcl9heDJbaV0gPC0gY3QuaSRlc3RpbWF0ZQ0KICBwdl9heDJbaV0gPC0gY3QuaSRwLnZhbHVlDQp9DQoNCmNvcl90YWIgPC0gZGF0YS5mcmFtZShlbnYgPSBuYW1lcyhlbnZfY29yciksDQogICAgICAgICAgICAgICAgICAgICAgb3JkX2F4MSA9IGNvcl9heDEsDQogICAgICAgICAgICAgICAgICAgICAgcHZhbF9heDEgPSBwdl9heDEsDQogICAgICAgICAgICAgICAgICAgICAgb3JkX2F4MiA9IGNvcl9heDIsDQogICAgICAgICAgICAgICAgICAgICAgcHZhbF9heDIgPSBwdl9heDIpDQoNCmNvcl90YWINCg0Kd3JpdGUuY3N2KHggPSBjb3JfdGFiLA0KICAgICAgICAgIGZpbGUgPSBmaWxlLnBhdGgoZGF0YVBhdGgsImlucHV0cy9Mb0RlbnNfY29yLXRhYmxlX3I2XzIwMHJlcF9NTEQtYmF0aHkuY3N2IiksDQogICAgICAgICAgcm93Lm5hbWVzID0gRkFMU0UpDQpgYGANCiMjIyBEb3QgY2hhcnQgdG8gY2hlY2sgZm9yIGdhcHMgaW4gY29ycmVsYXRpb24NCg0KYGBge3J9DQpjb3JfYTFfc29ydDwtY29yX3RhYiU+JQ0KICBtdXRhdGUoYWJzX29yZF9heDE9YWJzKG9yZF9heDEpLA0KICAgICAgICAgYWJzX29yZF9heDI9YWJzKG9yZF9heDIpKSAlPiUNCiAgYXJyYW5nZShkZXNjKGFic19vcmRfYXgxKSkNCg0KY29yX2EyX3NvcnQ8LWNvcl90YWIlPiUNCiAgbXV0YXRlKGFic19vcmRfYXgxPWFicyhvcmRfYXgxKSwNCiAgICAgICAgIGFic19vcmRfYXgyPWFicyhvcmRfYXgyKSkgJT4lDQogIGFycmFuZ2UoZGVzYyhhYnNfb3JkX2F4MikpDQoNCmRvdGNoYXJ0KGNvcl9hMV9zb3J0JGFic19vcmRfYXgxLCBtYWluPSJBYnNvbHV0ZSAoKy8tKSBjb3JyZWxhdGlvbnMgYmV0d2VlbiBlbnZWYXJzIGFuZCBnbm1kcyBheGlzIDEiKQ0KDQpjb3JfY3V0PC0wLjM3ICNkZWNpZGUNCg0KY29yX3NlbDwtc3Vic2V0KGNvcl9hMV9zb3J0LGFic19vcmRfYXgxPmNvcl9jdXQpDQpjb3Jfc2VsDQpgYGANCg0KDQojIyMgU2VsIGVudiB2YXIgKHRvcCBjb3JyKQ0KYGBge3J9DQplbnZfb3MgPC0gZW52WywgY29yX3NlbCRlbnZdDQplbnZfb3MNCnN0cihlbnZfb3MpDQoNCmBgYA0KDQoNCiMjIyBPcmRpc3VyZnMgdG9wIGNvcnINCmBgYHtyfQ0Kb3Jkc3JmcyA8LSBsaXN0KGxlbmd0aCA9IG5jb2woZW52X29zKSkNCg0KZm9yIChpIGluIHNlcShuY29sKGVudl9vcykpKSB7DQogIG9zLmkgPC0gZ2dfb3JkaXN1cmYob3JkID0gb3JkLA0KICAgICAgICAgICAgICAgICAgICAgIGVudi52YXIgPSBlbnZfb3NbLCBpXSwNCiAgICAgICAgICAgICAgICAgICAgICBwdC5zaXplID0gMSwNCiAgICAgICAgICAgICAgICAgICAgICAjIGJpbndpZHRoID0gMC4wNSwNCiAgICAgICAgICAgICAgICAgICAgICB2YXIubGFiZWwgPSBuYW1lcyhlbnZfb3MpW2ldLA0KICAgICAgICAgICAgICAgICAgICAgIGdlbi50ZXh0LnNpemUgPSAxMCwNCiAgICAgICAgICAgICAgICAgICAgICB0aXRsZS50ZXh0LnNpemUgPSAxNSwNCiAgICAgICAgICAgICAgICAgICAgICBsZWcudGV4dC5zaXplID0gMTApDQogIA0KICBvcmRzcmZzW1tpXV0gPC0gb3MuaSRwbG90DQp9DQoNCm9yZHNyZnNfcGx0IDwtIGdnYXJyYW5nZShwbG90bGlzdCA9IG9yZHNyZnMsDQogICAgICAgICAgICAgICAgICAgICAgICAgbnJvdyA9IDQsDQogICAgICAgICAgICAgICAgICAgICAgICAgbmNvbCA9IDMpDQoNCm9yZHNyZnNfcGx0DQpgYGANCiMjIyMjIFNhdmUgc29tZSBvdXRwdXRzDQpgYGB7cn0NCmdnZXhwb3J0KG9yZHNyZnNfcGx0LA0KICAgICAgICAgIGZpbGVuYW1lID0gZmlsZS5wYXRoKGRhdGFQYXRoLCJvdXRwdXRzL0xvRGVuc19vcmRpc3VyZnNfdG9wX2NvcnIucG5nIiksDQogICAgICAgICAgd2lkdGggPSAyMDAwLA0KICAgICAgICAgIGhlaWdodCA9IDI1MDApDQoNCmBgYA0KDQoNCiMjIyBTZWwgZW52IHZhciAobWFudWFsKQ0KYGBge3J9DQplbnZfb3NfbSA8LSBlbnZbLGMoIlRtZWFuX1JvYmluc29uIiwgI3RvcCBjb3JyDQogICAgICAgICAgICAgICAgICAic2FsdF9tYXgiLCAjdG9wIGNvcnINCiAgICAgICAgICAgICAgICAgICJTbWF4X1JvYmluc29uIiwgI2NvbXBhcmlzb24gdG8gdG9wIGNvcnINCiAgICAgICAgICAgICAgICAgICJzd0RlbnNSb2JfYXZzIiwgI3RvcCBjb3JyDQogICAgICAgICAgICAgICAgICAiQk8yMl9pY2Vjb3Zlcmx0bWF4X3NzIiwjdG9wIGNvcnIgYXgyDQogICAgICAgICAgICAgICAgICAiQk8yMl9pY2Vjb3Zlcm1lYW5fc3MiLCN0b3AgY29ycg0KICAgICAgICAgICAgICAgICAgIkJPMjJfZGlzc294bWVhbl9iZG1lYW4iLCN0b3AgY29ycg0KICAgICAgICAgICAgICAgICAgIyJCTzIyX2NhcmJvbnBoeXRvbHRtaW5fYmRtZWFuIiwjdG9wIGNvcnIgLSBubyBjbGVhciBncmFkaWVudCBpbiBvcmRpc3VyZg0KICAgICAgICAgICAgICAgICAgIkJPMjJfcHBsdG1pbl9zcyIsICN0b3AgY29ycg0KICAgICAgICAgICAgICAgICAgIlgueSIsICNjb21wYXJpc29uIHRvIFkNCiAgICAgICAgICAgICAgICAgICJZIiwgI3RvcCBjb3JyDQogICAgICAgICAgICAgICAgICAic3BkX3N0ZCIsICN0b3AgY29yciBheDIgKGJsZW5kZWQgbW9kZWwpDQogICAgICAgICAgICAgICAgICAiQ1NwZHNkX1JvYmluc29uIiwgI2NvbXBhcmlzb24gdG8gdG9wIGNvcnIgYXgyIChibGVuZGVkIG1vZGVsKQ0KICAgICAgICAgICAgICAgICAgIm11ZCIsICNoaWdoZXN0IHNlZCB2YXIgYXgxICsgY29ycg0KICAgICAgICAgICAgICAgICAgImdyYXZlbCIsI2hpZ2hlc3Qgc2VkIHZhciBheDEgLSBjb3JyDQogICAgICAgICAgICAgICAgICAiQk8yMl9zaWxpY2F0ZWx0bWF4X2JkbWVhbiIsICNqdXN0IHVuZGVyIHRvcCBjb3JyIGF4MQ0KICAgICAgICAgICAgICAgICAgImJhdGh5IiAjaW50dWl0aXZlIGZvciBjb21wYXJpc29ucw0KICAgICAgICAgICAgICAgICldDQplbnZfb3NfbQ0Kc3RyKGVudl9vc19tKQ0KDQpgYGANCg0KDQojIyMgT3JkaXN1cmZzIHRvcCBjb3JyDQpgYGB7cn0NCm9yZHNyZnNfbSA8LSBsaXN0KGxlbmd0aCA9IG5jb2woZW52X29zX20pKQ0KDQpmb3IgKGkgaW4gc2VxKG5jb2woZW52X29zX20pKSkgew0KICBvcy5pX20gPC0gZ2dfb3JkaXN1cmYob3JkID0gb3JkLA0KICAgICAgICAgICAgICAgICAgICAgIGVudi52YXIgPSBlbnZfb3NfbVssIGldLA0KICAgICAgICAgICAgICAgICAgICAgIHB0LnNpemUgPSAxLA0KICAgICAgICAgICAgICAgICAgICAgICMgYmlud2lkdGggPSAwLjA1LA0KICAgICAgICAgICAgICAgICAgICAgIHZhci5sYWJlbCA9IG5hbWVzKGVudl9vc19tKVtpXSwNCiAgICAgICAgICAgICAgICAgICAgICBnZW4udGV4dC5zaXplID0gMTAsDQogICAgICAgICAgICAgICAgICAgICAgdGl0bGUudGV4dC5zaXplID0gMTUsDQogICAgICAgICAgICAgICAgICAgICAgbGVnLnRleHQuc2l6ZSA9IDEwKQ0KICANCiAgb3Jkc3Jmc19tW1tpXV0gPC0gb3MuaV9tJHBsb3QNCn0NCg0Kb3Jkc3Jmc19wbHRfbSA8LSBnZ2FycmFuZ2UocGxvdGxpc3QgPSBvcmRzcmZzX20sDQogICAgICAgICAgICAgICAgICAgICAgICAgbnJvdyA9IDQsDQogICAgICAgICAgICAgICAgICAgICAgICAgbmNvbCA9IDQpDQoNCm9yZHNyZnNfcGx0X20NCmBgYA0KDQojIyMjIyBTYXZlIHNvbWUgb3V0cHV0cw0KYGBge3J9DQpnZ2V4cG9ydChvcmRzcmZzX3BsdF9tLA0KICAgICAgICAgIGZpbGVuYW1lID0gZmlsZS5wYXRoKGRhdGFQYXRoLCJvdXRwdXRzL0xvRGVuc19vcmRpc3VyZnNfbWFuX3NlbC5wbmciKSwNCiAgICAgICAgICB3aWR0aCA9IDIwMDAsDQogICAgICAgICAgaGVpZ2h0ID0gMjAwMCkNCg0KYGBgDQoNCg0KDQoNCiMjIyBFbnZmaXQNCmBgYHtyfQ0KIyMgU2VsZWN0IGlmIGFueSB2YXIgc2hvdWxkIGJlIGV4Y2x1ZGVkIGZyb20gZW52Zml0IChtYWtlcyBsZXNzIGJ1c3kgdG8gcmVhZCkNCmVudl9vc19tX2VudmZpdDwtZW52X29zX20gWyxjKCJUbWVhbl9Sb2JpbnNvbiIsICN0b3AgY29ycg0KICAgICAgICAgICAgICAgICAgInNhbHRfbWF4IiwgI3RvcCBjb3JyDQogICAgICAgICAgICAgICAgICAiU21heF9Sb2JpbnNvbiIsICNjb21wYXJpc29uIHRvIHRvcCBjb3JyDQogICAgICAgICAgICAgICAgICAic3dEZW5zUm9iX2F2cyIsICN0b3AgY29ycg0KICAgICAgICAgICAgICAgICAgIkJPMjJfaWNlY292ZXJsdG1heF9zcyIsI3RvcCBjb3JyIGF4Mg0KICAgICAgICAgICAgICAgICAgIkJPMjJfaWNlY292ZXJtZWFuX3NzIiwjdG9wIGNvcnINCiAgICAgICAgICAgICAgICAgICJCTzIyX2Rpc3NveG1lYW5fYmRtZWFuIiwjdG9wIGNvcnINCiAgICAgICAgICAgICAgICAgICMiQk8yMl9jYXJib25waHl0b2x0bWluX2JkbWVhbiIsI3RvcCBjb3JyIC0gbm8gY2xlYXIgZ3JhZGllbnQgaW4gb3JkaXN1cmYNCiAgICAgICAgICAgICAgICAgICJCTzIyX3BwbHRtaW5fc3MiLCAjdG9wIGNvcnINCiAgICAgICAgICAgICAgICAgICJYLnkiLCAjY29tcGFyaXNvbiB0byBZDQogICAgICAgICAgICAgICAgICAiWSIsICN0b3AgY29ycg0KICAgICAgICAgICAgICAgICAjICJzcGRfc3RkIiwgI3RvcCBjb3JyIGF4MiAoYmxlbmRlZCBtb2RlbCkNCiAgICAgICAgICAgICAgICAgICJDU3Bkc2RfUm9iaW5zb24iLCAjY29tcGFyaXNvbiB0byB0b3AgY29yciBheDIgKGJsZW5kZWQgbW9kZWwpDQogICAgICAgICAgICAgICAgICAibXVkIiwgI2hpZ2hlc3Qgc2VkIHZhciBheDEgKyBjb3JyDQogICAgICAgICAgICAgICAgICAiZ3JhdmVsIiwjaGlnaGVzdCBzZWQgdmFyIGF4MSAtIGNvcnINCiAgICAgICAgICAgICAgICAgICJCTzIyX3NpbGljYXRlbHRtYXhfYmRtZWFuIiwgI2p1c3QgdW5kZXIgdG9wIGNvcnIgYXgxDQogICAgICAgICAgICAgICAgICAiYmF0aHkiICNpbnR1aXRpdmUgZm9yIGNvbXBhcmlzb25zDQogICAgICAgICAgICAgICAgDQogIA0KKV0NCg0KY29sbmFtZXMoZW52X29zX21fZW52Zml0KTwtYygiVCIsICN0b3AgY29ycg0KICAgICAgICAgICAgICAgICAgInNNeF9fX18iLCAjdG9wIGNvcnINCiAgICAgICAgICAgICAgICAgICJTbWF4UiIsICNjb21wYXJpc29uIHRvIHRvcCBjb3JyDQogICAgICAgICAgICAgICAgICAic3dEZW5zUiIsICN0b3AgY29ycg0KICAgICAgICAgICAgICAgICAgImljZWNvdm1heCIsI3RvcCBjb3JyIGF4Mg0KICAgICAgICAgICAgICAgICAgImljZWNvdmF2IiwjdG9wIGNvcnINCiAgICAgICAgICAgICAgICAgICJkaXNzb3hhdiIsI3RvcCBjb3JyDQogICAgICAgICAgICAgICAgICAjIkJPMjJfY2FyYm9ucGh5dG9sdG1pbl9iZG1lYW4iLCN0b3AgY29yciAtIG5vIGNsZWFyIGdyYWRpZW50IGluIG9yZGlzdXJmDQogICAgICAgICAgICAgICAgICAicHBsdG1pbiIsICN0b3AgY29ycg0KICAgICAgICAgICAgICAgICAgIlgiLCAjY29tcGFyaXNvbiB0byBZDQogICAgICAgICAgICAgICAgICAiWSIsICN0b3AgY29ycg0KICAgICAgICAgICAgICAgICAjICJzcGRfc3RkIiwgI3RvcCBjb3JyIGF4MiAoYmxlbmRlZCBtb2RlbCkNCiAgICAgICAgICAgICAgICAgICJDU3NkIiwgI2NvbXBhcmlzb24gdG8gdG9wIGNvcnIgYXgyIChibGVuZGVkIG1vZGVsKQ0KICAgICAgICAgICAgICAgICAgIm11ZCIsICNoaWdoZXN0IHNlZCB2YXIgYXgxICsgY29ycg0KICAgICAgICAgICAgICAgICAgImdyYXZlbCIsI2hpZ2hlc3Qgc2VkIHZhciBheDEgLSBjb3JyDQogICAgICAgICAgICAgICAgICAiU2lMdG1heCIsICNqdXN0IHVuZGVyIHRvcCBjb3JyIGF4MQ0KICAgICAgICAgICAgICAgICAgImJhdGh5IiAjaW50dWl0aXZlIGZvciBjb21wYXJpc29ucw0KICAgICAgICAgICAgICAgICApDQoNCiMjIEVudmZvdCBwbG90DQpnZ19lbnZmaXQob3JkID0gb3JkLA0KICAgICAgICAgIGVudiA9IGVudl9vc19tX2VudmZpdCwNCiAgICAgICAgICBwdC5zaXplID0gMSkNCg0KIyMgRW52Zml0IGFuYWx5c2lzDQoNCmVmIDwtIGVudmZpdChvcmQgPSBvcmQsDQogICAgICAgICAgICAgZW52ID0gZW52X29zX21fZW52Zml0LA0KICAgICAgICAgICAgICMgbmEucm0gPSBUUlVFDQogICAgICAgICAgICAgKQ0KZWZERiA8LSBhcy5kYXRhLmZyYW1lKHNjb3JlcyhlZiwNCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgZGlzcGxheSA9ICJ2ZWN0b3JzIikpDQpgYGANCg0KIyMjIyMgU2F2ZSB0aGUgcGxvdA0KYGBge3J9DQpnZ3NhdmUoZmlsZW5hbWUgPSBmaWxlLnBhdGgoZGF0YVBhdGgsIm91dHB1dHMvTG9EZW5zX0VudkZpdF9tYW5fc2VsX2Nsbi5wbmciKSwNCiAgICAgICBkZXZpY2UgPSAicG5nIiwNCiAgICAgICBkcGk9MzAwICkNCmBgYA0KIyMjIENhdGVnb3JpY2FsIGVudlZhciB2aXN1YWxpc2VkIG9uIHRoZSBtZHNwbG90cw0KDQojIyMjIHVzZSBkYXRhc2V0IHdpdGggY2F0ZWdvcmljYWwgdmFyIGluY2x1ZGVkDQoNCmBgYHtyfQ0KZW52X3ZpczwtZW52DQplbnZfdmlzJGdubWRzMSA8LSBvdHVfNiRnbm1kczENCmVudl92aXMkZ25tZHMyIDwtIG90dV82JGdubWRzMg0KYGBgDQoNCiMjIyMgZ25tZHMgdyBtbGQgbWVhbiAtIGJhdGh5DQpgYGB7cn0NCnBfbWxkIDwtIGdncGxvdChkYXRhID0gZW52X3ZpcywNCiAgICAgICAgICAgICAgICAgICAgIGFlcyh4ID0gZ25tZHMxLA0KICAgICAgICAgICAgICAgICAgICAgICAgIHkgPSBnbm1kczIpKSArDQogIHRoZW1lX2NsYXNzaWMoKSArDQogIGNvb3JkX2ZpeGVkKCkgKw0KICBnZ3RpdGxlKCJHTk1EUyBjb2xvdXJlZCBieSBwcm94aW1pdHkgdG8gbWl4ZWQgbGF5ZXIgZGVwdGgiLA0KICAgICAgICAgIHN1YnRpdGxlID0gIkZpcnN0IHJ1biIpICsNCiAgZ2VvbV9wb2ludChhZXMoY29sb3VyID0gTUxEbWVhbl9iYXRoeSkpICsNCiAgZ2VvbV92bGluZSh4aW50ZXJjZXB0ID0gMCwNCiAgICAgICAgICAgICBsaW5ldHlwZSA9IDIsDQogICAgICAgICAgICAgY29sb3VyID0gImxpZ2h0Z3JheSIpICsNCiAgZ2VvbV9obGluZSh5aW50ZXJjZXB0ID0gMCwNCiAgICAgICAgICAgICBsaW5ldHlwZSA9IDIsDQogICAgICAgICAgICAgY29sb3VyID0gImxpZ2h0Z3JheSIpDQoNCnBfbWxkDQpgYGANCiMjIyMjIFNhdmUgdGhlIHBsb3QNCmBgYHtyfQ0KZ2dzYXZlKGZpbGVuYW1lID0gZmlsZS5wYXRoKGRhdGFQYXRoLCJvdXRwdXRzL0xvRGVuc19nbm1kc19NTERtZWFuQmF0aHkucG5nIiksDQogICAgICAgZGV2aWNlID0gInBuZyIsDQogICAgICAgZHBpPTMwMCApDQpgYGANCiMjIyMgZ25tZHMgdyBzZWRjbGFzcw0KIyMjIyMgZGVhbCB3aXRoIHNlZGNsYXNzIGNvZGVzDQpgYGB7cn0NCmVudl92aXM8LSBlbnZfdmlzICU+JSANCiAgbXV0YXRlKA0KICAgIHNlZGNsYXNzTmFtZSA9IGNhc2Vfd2hlbigNCiAgICAgICAgICAgIHNlZGNsYXNzID09ICIxIiB+ICJTZWRDb3ZlclIiLA0KICAgICAgICAgICAgc2VkY2xhc3MgPT0gIjUiIH4gIlJvY2siLA0KICAgICAgICAgICAgc2VkY2xhc3MgPT0gIjIwIiB+ICJNdWQiLA0KICAgICAgICAgICAgc2VkY2xhc3MgPT0gIjIxIiB+ICJNd0Jsb2NrIiwNCiAgICAgICAgICAgIHNlZGNsYXNzID09ICI0MCIgfiAic011ZCIsDQogICAgICAgICAgICBzZWRjbGFzcyA9PSAiODAiIH4gIm1TYW5kIiwNCiAgICAgICAgICAgIHNlZGNsYXNzID09ICIxMDAiIH4gIlNhbmQiLA0KICAgICAgICAgICAgc2VkY2xhc3MgPT0gIjExMCIgfiAiZ011ZCIsDQogICAgICAgICAgICBzZWRjbGFzcyA9PSAiMTE1IiB+ICJnc011ZCIsDQogICAgICAgICAgICBzZWRjbGFzcyA9PSAiMTIwIiB+ICJnbVNhbmQiLA0KICAgICAgICAgICAgc2VkY2xhc3MgPT0gIjEzMCIgfiAiZ1NhbmQiLA0KICAgICAgICAgICAgc2VkY2xhc3MgPT0gIjE1MCIgfiAiTVNHIiwNCiAgICAgICAgICAgIHNlZGNsYXNzID09ICIxNjAiIH4gInNHcmF2ZWwiLA0KICAgICAgICAgICAgc2VkY2xhc3MgPT0gIjE3MCIgfiAiR3JhdmVsIiwNCiAgICAgICAgICAgIHNlZGNsYXNzID09ICIxNzUiIH4gIkdyYXZCbG9jayIsDQogICAgICAgICAgICBzZWRjbGFzcyA9PSAiMTg1IiB+ICJTR0JtaXgiLA0KICAgICAgICAgICAgc2VkY2xhc3MgPT0gIjIwNSIgfiAiUy9Nd0IiLA0KICAgICAgICAgICAgc2VkY2xhc3MgPT0gIjIwNiIgfiAiUy9Nd0cvQiIsDQogICAgICAgICAgICBzZWRjbGFzcyA9PSAiMjE1IiB+ICJTR0JhbHQiLA0KICAgICAgICAgICAgc2VkY2xhc3MgPT0gIjMwMCIgfiAiSGFyZFNlZCIsDQogICAgICAgICAgICBzZWRjbGFzcyA9PSAiNTAwIiB+ICJCaW9nZW5pYyINCiAgICAgICAgICAgICAgICAgICAgICAgIA0KICAgICkNCiAgKQ0KYGBgDQoNCiMjIyMjIyBjb2xvdXIgcGFsZXR0ZSB0byBjb3BlIHdpdGggdXAgdG8gMjUgY2F0ZWdvcmljYWwgY29sb3Vycw0KYGBge3J9DQpjMjUgPC0gYygNCiAgImRvZGdlcmJsdWUyIiwgIiNFMzFBMUMiLCAjIHJlZA0KICAiZ3JlZW40IiwNCiAgIiM2QTNEOUEiLCAjIHB1cnBsZQ0KICAiI0ZGN0YwMCIsICMgb3JhbmdlDQogICJibGFjayIsICJnb2xkMSIsDQogICJza3libHVlMiIsICIjRkI5QTk5IiwgIyBsdCBwaW5rDQogICJwYWxlZ3JlZW4yIiwNCiAgIiNDQUIyRDYiLCAjIGx0IHB1cnBsZQ0KICAiI0ZEQkY2RiIsICMgbHQgb3JhbmdlDQogICJncmF5NzAiLCAia2hha2kyIiwNCiAgIm1hcm9vbiIsICJvcmNoaWQxIiwgImRlZXBwaW5rMSIsICJibHVlMSIsICJzdGVlbGJsdWU0IiwNCiAgImRhcmt0dXJxdW9pc2UiLCAiZ3JlZW4xIiwgInllbGxvdzQiLCAieWVsbG93MyIsDQogICJkYXJrb3JhbmdlNCIsICJicm93biINCikNCmBgYA0KDQoNCg0KDQpgYGB7cn0NCg0KcF9zZWQgPC0gZ2dwbG90KGRhdGEgPSBlbnZfdmlzLA0KICAgICAgICAgICAgICAgICAgICAgYWVzKHggPSBnbm1kczEsDQogICAgICAgICAgICAgICAgICAgICAgICAgeSA9IGdubWRzMikpICsNCiAgdGhlbWVfY2xhc3NpYygpICsNCiAgY29vcmRfZml4ZWQoKSArDQogIGdndGl0bGUoIkdOTURTIGNvbG91cmVkIGJ5IHNlZGltZW50IGNsYXNzIiwNCiAgICAgICAgICBzdWJ0aXRsZSA9ICJGaXJzdCBydW4iKSArDQogIGdlb21fcG9pbnQoYWVzKGNvbG91ciA9IGZhY3RvcihzZWRjbGFzc05hbWUpKSkgKw0KICAgc2NhbGVfY29sb3VyX21hbnVhbCh2YWx1ZXM9YzI1KSsNCiAgZ2VvbV92bGluZSh4aW50ZXJjZXB0ID0gMCwNCiAgICAgICAgICAgICBsaW5ldHlwZSA9IDIsDQogICAgICAgICAgICAgY29sb3VyID0gImxpZ2h0Z3JheSIpICsNCiAgZ2VvbV9obGluZSh5aW50ZXJjZXB0ID0gMCwNCiAgICAgICAgICAgICBsaW5ldHlwZSA9IDIsDQogICAgICAgICAgICAgY29sb3VyID0gImxpZ2h0Z3JheSIpKw0KICBndWlkZXMoY29sb3VyPWd1aWRlX2xlZ2VuZChuY29sPTIpKQ0KDQpwX3NlZA0KYGBgDQojIyMjIyBTYXZlIHRoZSBwbG90DQpgYGB7cn0NCmdnc2F2ZShmaWxlbmFtZSA9IGZpbGUucGF0aChkYXRhUGF0aCwib3V0cHV0cy9Mb0RlbnNfZ25tZHNfc2VkY2xhc3MucG5nIiksDQogICAgICAgZGV2aWNlID0gInBuZyIsDQogICAgICAgZHBpPTMwMCApDQpgYGANCiMjIyMgZ25tZHMgdyBsYW5zY2FwZQ0KIyMjIyMgZGVhbCB3aXRoIHNlZGNsYXNzIGNvZGVzDQpgYGB7cn0NCmVudl92aXM8LSBlbnZfdmlzICU+JSANCiAgbXV0YXRlKA0KICAgIGxhbmRzY2FwZU5hbWUgPSBjYXNlX3doZW4oDQogICAgICAgICAgICBsYW5kc2NhcGUgPT0gIjEiIH4gIlN0cmFuZGZsYXQiLA0KICAgICAgICAgICAgbGFuZHNjYXBlID09ICIyMSIgfiAiQ29udFNsb3BlIiwNCiAgICAgICAgICAgIGxhbmRzY2FwZSA9PSAiMjIiIH4gIkNhbnlvbiIsDQogICAgICAgICAgICBsYW5kc2NhcGUgPT0gIjMxIiB+ICJWYWxsZXkiLA0KICAgICAgICAgICAgbGFuZHNjYXBlID09ICIzMiIgfiAiRmpvcmQiLA0KICAgICAgICAgICAgbGFuZHNjYXBlID09ICI0MSIgfiAiRGVlcFBsYWluIiwNCiAgICAgICAgICAgIGxhbmRzY2FwZSA9PSAiNDIiIH4gIlNsb3BlUGxhaW4iLA0KICAgICAgICAgICAgbGFuZHNjYXBlID09ICI0MyIgfiAiU2hlbGZQbGFpbiIsDQogICAgICAgICAgICBsYW5kc2NhcGUgPT0gIjQzMSIgfiAic2hhbGxvd1ZhbGxleSINCiAgICApDQogICkNCmBgYA0KDQoNCmBgYHtyfQ0KDQpwX2xhbmQgPC0gZ2dwbG90KGRhdGEgPSBlbnZfdmlzLA0KICAgICAgICAgICAgICAgICAgICAgYWVzKHggPSBnbm1kczEsDQogICAgICAgICAgICAgICAgICAgICAgICAgeSA9IGdubWRzMikpICsNCiAgdGhlbWVfY2xhc3NpYygpICsNCiAgY29vcmRfZml4ZWQoKSArDQogIGdndGl0bGUoIkdOTURTIGNvbG91cmVkIGJ5IGxhbmRzY2FwZSBjbGFzcyIsDQogICAgICAgICAgc3VidGl0bGUgPSAiRmlyc3QgcnVuIikgKw0KICBnZW9tX3BvaW50KGFlcyhjb2xvdXIgPSBmYWN0b3IobGFuZHNjYXBlTmFtZSkpKSArDQogIGdlb21fdmxpbmUoeGludGVyY2VwdCA9IDAsDQogICAgICAgICAgICAgbGluZXR5cGUgPSAyLA0KICAgICAgICAgICAgIGNvbG91ciA9ICJsaWdodGdyYXkiKSArDQogIGdlb21faGxpbmUoeWludGVyY2VwdCA9IDAsDQogICAgICAgICAgICAgbGluZXR5cGUgPSAyLA0KICAgICAgICAgICAgIGNvbG91ciA9ICJsaWdodGdyYXkiKSsNCiAgZ3VpZGVzKGNvbG91cj1ndWlkZV9sZWdlbmQobmNvbD0yKSkNCg0KcF9sYW5kDQpgYGANCiMjIyMjIFNhdmUgdGhlIHBsb3QNCmBgYHtyfQ0KZ2dzYXZlKGZpbGVuYW1lID0gZmlsZS5wYXRoKGRhdGFQYXRoLCJvdXRwdXRzL0xvRGVuc19nbm1kc19sYW5kc2NhcGUucG5nIiksDQogICAgICAgZGV2aWNlID0gInBuZyIsDQogICAgICAgZHBpPTMwMCApDQpgYGANCg0KIyMjIyBHbm1kcyB3IGdtb3JwaA0KYGBge3J9DQoNCnBfZ21vIDwtIGdncGxvdChkYXRhID0gZW52X3ZpcywNCiAgICAgICAgICAgICAgICAgICAgIGFlcyh4ID0gZ25tZHMxLA0KICAgICAgICAgICAgICAgICAgICAgICAgIHkgPSBnbm1kczIpKSArDQogIHRoZW1lX2NsYXNzaWMoKSArDQogIGNvb3JkX2ZpeGVkKCkgKw0KICBnZ3RpdGxlKCJHTk1EUyBjb2xvdXJlZCBieSBsYW5kc2NhcGUgY2xhc3MiLA0KICAgICAgICAgIHN1YnRpdGxlID0gIkZpcnN0IHJ1biIpICsNCiAgZ2VvbV9wb2ludChhZXMoY29sb3VyID0gZmFjdG9yKGdtb3JwaCkpKSArDQogIGdlb21fdmxpbmUoeGludGVyY2VwdCA9IDAsDQogICAgICAgICAgICAgbGluZXR5cGUgPSAyLA0KICAgICAgICAgICAgIGNvbG91ciA9ICJsaWdodGdyYXkiKSArDQogIGdlb21faGxpbmUoeWludGVyY2VwdCA9IDAsDQogICAgICAgICAgICAgbGluZXR5cGUgPSAyLA0KICAgICAgICAgICAgIGNvbG91ciA9ICJsaWdodGdyYXkiKSsNCiAgZ3VpZGVzKGNvbG91cj1ndWlkZV9sZWdlbmQobmNvbD0yKSkNCg0KcF9nbW8NCmBgYA0KIyMjIyMgU2F2ZSB0aGUgcGxvdA0KYGBge3J9DQpnZ3NhdmUoZmlsZW5hbWUgPSBmaWxlLnBhdGgoZGF0YVBhdGgsIm91dHB1dHMvTG9EZW5zX2dubWRzX2dtb3JwaC5wbmciKSwNCiAgICAgICBkZXZpY2UgPSAicG5nIiwNCiAgICAgICBkcGk9MzAwICkNCmBgYA0KDQpgYGB7cn0NCmNhdF92YXJfcGxvdHM8LXBfbWxkK3Bfc2VkK3BfbGFuZCtwX2dtbw0KYGBgDQojIyMjIyBTYXZlIHRoZSBwbG90DQpgYGB7cn0NCmdnZXhwb3J0KGNhdF92YXJfcGxvdHMsDQogICAgICAgICAgZmlsZW5hbWUgPSBmaWxlLnBhdGgoZGF0YVBhdGgsIm91dHB1dHMvTG9EZW5zX2dubWRzX2NhdHZhci5wbmciKSwNCiAgICAgICAgICB3aWR0aCA9IDEwMDAsDQogICAgICAgICAgaGVpZ2h0ID0gODAwKQ0KYGBgDQoNCiMjIyMgU2FtcGxlIGlkZW50aWZpY2F0aW9uIGluIHRoZSBtZHNwbG90DQoNCmBgYHtyfQ0KcF9nbW8gPC0gZ2dwbG90KGRhdGEgPSBlbnZfdmlzLA0KICAgICAgICAgICAgICAgICAgICAgYWVzKHggPSBnbm1kczEsDQogICAgICAgICAgICAgICAgICAgICAgICAgeSA9IGdubWRzMikpICsNCiAgdGhlbWVfY2xhc3NpYygpICsNCiAgY29vcmRfZml4ZWQoKSArDQogIGdndGl0bGUoIkdOTURTIGNvbG91cmVkIGJ5IGxhbmRzY2FwZSBjbGFzcyIsDQogICAgICAgICAgc3VidGl0bGUgPSAiRmlyc3QgcnVuIikgKw0KICBnZW9tX3BvaW50KGFlcyhjb2xvdXIgPSBmYWN0b3IoU2FtcElEKSkpICsNCiAgZ2VvbV92bGluZSh4aW50ZXJjZXB0ID0gMCwNCiAgICAgICAgICAgICBsaW5ldHlwZSA9IDIsDQogICAgICAgICAgICAgY29sb3VyID0gImxpZ2h0Z3JheSIpICsNCiAgZ2VvbV9obGluZSh5aW50ZXJjZXB0ID0gMCwNCiAgICAgICAgICAgICBsaW5ldHlwZSA9IDIsDQogICAgICAgICAgICAgY29sb3VyID0gImxpZ2h0Z3JheSIpKw0KICBndWlkZXMoY29sb3VyPSJub25lIiwgc2l6ZT0ibm9uZSIpDQoNCmdncGxvdGx5KHBfZ21vKQ0KYGBgDQoNCg==